AD-A183  892 


T 


AFVAL-TR-87- 3034 


1.ASER  INDUCED  EMP  AT  10.6  MICRONS 


(Hen  Dahlbacka  and  John  Guillory 

Plasma  Research  Corp. 

2900  Main  Street 
Alameda,  CA  94501 


May  1987 


Final  Report  for  Period  January  1984  -  June  1986 


Approved  for  public  release;  distribution  unlimited. 


FLIGHT  DYNAMICS  LABORATORY 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 

AIR  FORCE  SYSTEMS  COMMAND 

WRIGHT- PATTERSON  AIR  FORCE  BASE,  OHIO  45433-6553 


87  8  18  055 


u'<r.aia<F.wn:om:K.MnTM7.Tcn 


REPORT  DOCUMENTATION  PAGE 


2 a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 
PRC  FR  86-01 


6#  NAME  OF  PERFORMING  ORGANIZATION 

Plasma  Research  Corp. 

6b  OFFICE  SYMBOL 
(If  applicable) 

6c  ADDRESS  (Oty,  Stite,  tnd  ZIP  Code) 

2900  Main  St. 

Alameda,  CA  94501 

8«.  NAME  OF  FUNDING /SPONSORING 

8b  OFFICE  SYMBOL 

ORGANIZATION 

(If  applicable) 

AFWAL  FDL  AFSC 

FIESL 

8c.  ADDRESS  (City.  Stite,  tnd  ZIP  Code) 

Wright-Patterson  AFB,  Ohio 
45433-6553 

lb.  RESTRICTIVE  MARKINGS 
N/A 


3  DISTRIBUTION /AVAILABILITY  OF  REPORT 
Approved  for  public  release, 
dlscributlon  unlimited. 


S  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 
AFWAL-TR-87-3034 


7a  NAME  OF  MONITORING  ORGANIZATION 
Flight  Dynamics  Laboratory  (AFWAL/FI) 


7b  ADDRESS  (C/ty,  State,  and  ZIP  Cod*) 

Air  Force  Wright  Aeronautical  Laboratories 
Wright-Patterson  AFB,  OH  45433-6553 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F33615-85-C-3405 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 

62201F 


PROJECT 

TASK 

NO. 

NO 

2402 

02 

11  TITLE  (Include  Security  Clissifkition) 

Laser  Induced  EMP  at  10.6  Microns  Unclassified 


12  PERSONAL  AUTHOR(S) 

Dahlbacka,  Glen  PhD;  Guillory,  John  PhD 


14  DATE  OF  REPORT  (Year,  Month,  Day) 

1987  Ma 


16  SUPPLEMENTARY  NOTATION 


‘  ~  c  'J 


FIELD 

GROUP 

20 

09 

COSATI  CODES 


SUB-GROUP 


18  SUBJECT  TERMS  (Continue  on  reverse  if  necessity  and  identify  by  block  number ) 

lLaser  Plasma,'  Laser^  Electromagnetic  Pulsed  EMP, 
'Scaling  Laws# ^Experiment  &  Theory,  Physics 


19  ABSTRACT  (Continue  on  reverse  if  necessity  *nd  identify  by  block  number) 

^>An  extensive  experimental  session  was  analyzed  (1200  shots)  to  investigate  the  behavior 
of  Laser  Induced  ElectroMagnetic  Pulse  (LIEMP)  parametrically  with  C0(Z)  laser  intensity 
within  a  decade  of  10  GW/^ta(2)  and  with  air  pressure  from  1  microtorr  to  a  torr. 

Measured  electric  fields  on  a  ground  plane  and  current  measured  3  cm  from  the  target  were 
compared  to  theory.  One  hundred  eV  suprathermal  electrons  (1.6  keV  maximum  energy)  create 
the  LIEMP  after  laser  light,  self-focussed  a  few  fold,  is  resonantly  absorbed.  The  model 
predicts  electric  field  scaling  as  Intensity (1/2)  as  observed,  but  current  scaling  is  not 
reproduced.  The  electric  fields  are  no  threat  for  C0(2)  and  shorter  wavelengths,  but  may 
pose  a  threat  for  millimeter  wayes.  Pressure  dependencies  are  reproduced  with  classical 
energy  loss  in  gasses.  . 


2b  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 

duNCLASSIFIEP/UNLIMITEO  □  SAME  AS  RPT  □  OTIC  USERS 


22#  NAME  OF  RESPONSIBLE  INDIVIDUAL 
Capt  Harold  D.  Burket 


DO  FORM  1473,  84  MAR  83  APR  edition  may  be  used  until  exhausted 

All  other  editions  are  obsolete 


21  ABSTRACT  SECURITY  CLASSIFICATION 
CLASSIFIED 


22b  TELEPHONE  (Include  Area  Code)  22c  OFFICE  SYMBOL 
(513)  257-7718  AFWAL/FIESL 


SECURITY  CLASSIFICATION  OF  ~HIS  PAGE 
UNCLASSIFIED 


TABLE  OF  CONTENTS 


Chapter  1  TRACE . . . 

Chapter  2  INTRODUCTION . . . 

Chapter  3  THEORETICAL  ANALYSIS . 

3.0  Introduction  to  the  Theoretical  Analysis... 

3.1  Thresholds  for  Various  Processes  Known  to 
Produce  Fast  Electrons  at  Laser  Targets.... 

3.2  Fllamentatlon:  Mechanisms,  Thresholds, 
Amplitudes,  and  Effects.................... 

3.3  Space-Charge  Layer  Near  an  Electron- 

Emitting  Target  in  Vacuum . 

3.4  Current  Density  of  Escaping  Fast  Electrons. 

3.5  Fast  Electron  Transit  Times  from  Isolated 

vs  Grounded  Targets . 

3.6  Radial  Spreading  of  the  Fast-Electron  Plume 

3.7  Fast  Electron  Slowing  Down . . 

3.8  Plume  Electro-Magnetic  Fields  and  the 

Motion  of  Secondary  Electrons  in  the 
Collisionless  Limit.... . 

Chapter  4  DATA  ANALYSIS..... . 

4.1  General  Waveform  Analysis . 

4.2  Data  Tabulation . . . 

4.3  Analysis  of  Tabulation-Intensity  Scaling... 

4.4  Analysis  of  Tabulation-Pressure  Scaling.... 

4.5  Experimental  Scaling  Summary.. . 

Chapter  5  SCALING  OF  OBSERVABLES  WITH  LASER 

INTENSITY,  WAVELENGTH,  AND  BACKGROUND. 

5.1  Intensity  Scaling................ . 


I 

4 

8 

8 

12 

19 

21 

33 

40 

50 

52 

54 

68 

68 

71 

76 

78 

80 


99  _ 

99  — 

C 


iii 


5*2  Wavelength  Scaling 
5.3  Pressure  Scaling*. 


Chapter  6  ACKNOWLEDGEMENTS 


Appendix  A  Hot-Electron  Transit  Times  in  a  Self- 

Consistent  Grounded-Dipole  Potential.... 

Appendix  B  Hot-Electron  Transit  Times  in  a  Monopole 
Space-Charge  Potential  Well... . 

Appendix  C  Effect  of  Local  Magnetic  Fields  on 

Electron  Plume  Motion . . . 

References. . . . 

Conversion  Chart . 


Chapter  1 


TRACE 

An  extensive  experimental  session  was  conducted  to  Investigate  the 
behavior  of  Laser  Induced  ElectroMagnetlc  Pulses  (LIEMP)  parametrically 
with  CO2  laser  intensity  within  a  decade  of  10  GW/cm^»  Parametric 
variations  were  also  made  with  background  air  pressure  from  1  microtorr  to 
a  torr.  Twelve  hundred  (1200)  data  shots  were  taken  with  measurements  of 
electric  fields  on  a  ground  plane  outside  a  glass  vacuum  chamber,  Faraday 
cup  (PCI)  measurements  with  the  vacuum  chamber  30  cm  from  the  Copper 
target  and  30  degrees  off  the  Incident  laser  vector,  and  current 
measurements  In  the  target  plane  were  made  with  a  Rogowski  style  monitor 
(FMM)  within  a  3  cm  diameter  of  the  target. 

Theoretical  work  leads  to  the  following  model  for  the  behavior. 

The  Incident  laser  light  Is  self-focussed  such  that  its  Intensity  Is 
Increased  by  a  few  fold  before  It  reached  the  critical  surface.  At  that 
point,  resonance  absorption  places  laser  energy  Into  the  electrons  that 
have  a  suprathermal  distribution  with  a  hot  temperature  of  200  eV  or  so 
depending  on  intensity.  Some  of  these  electrons  then  leave  the  target  as 
an  electron  beam,  creating  the  LIEMP  effect  as  well  as  observable  currents 
in  the  ground  plane  FMM  detector  and  a  FC1  detector.  The  most  energetic 
electrons  observed  were  1.6  keV,  with  more  typical  energies  of  600  eV 
observed  30  cm  from  the  plasma. 

Analysis  of  the  experimental  data  leads  to  well  founded  scaling 
rules  for  the  currents  at  le-6  torr  and  for  the  Faraday  cup  currents  at 


both  pressures 


3/2 

PMM-current  •  Intensity 

3/4 

FC1  current  «  Intensity 

FC1  energy  *  Intensity** 

(~  1/2) 

FC1  current  «  FMM  current  ' 

For  the  electric  fields  there  is  conflicting  data  with  respect  to 
the  theory.  At  le-6  torr  the  bulk  of  the  data  scales  as 

Field  ~  Intensity3^2  , 

which  is  not  easily  made  consistent  with  theory  (~  1/2);  however,  one  data 
point  at  the  very  lowest  intensity  could  be  construed  to  bring  the  scaling 
to  the  1/2  power,  consistent  with  the  3e-3  torr  data  which  experimentally 
scales  as 


Field  ~  Intensity*^2  . 

The  electric  field  observations  were  made  difficult  by  high 
frequency  components  at  low  pressure. 

The  scaling  as  1^^  Is  consistent  with  radiation  from  a  virtual 
cathode  oscillator,  but  theory  shows  that  these  high  frequency 
oscillations  are  probably  smaller  than  the  quasi-static  plume  fields  that 


radiate  on  the  order  of  20  times  more  in  these  experiments.  We  use  the 

quasi-static  plume  E-field  scaling  in  scaling  to  higher  laser  energies. 

1/3 

The  quasi-static  plume  field  scales  as  (wavelength)  at  a  constant 
intensity.  Thus  there  is  a  weak  supression  of  these  fields  at  shorter  wave¬ 
lengths.  This  means  that  electric  fields  created  by  CO^  lasers  and  shorter 
wavelength  lasers  have  little  system  effect  from  the  field.  Millimeter  wave 
lasers,  however  can  create  Megavolt  potentials  for  isolated  targets  in  space 
and  may  be  a  threat.  The  laser  must  create  a  plasma  for  these  effects  to 
become  apparent. 

The  scaling  of  field  with  pressure  Is  well  explained  by  classical  dE/dx 
of  the  electrons  in  the  background  gas  leading  to  reduced  charge  separations. 
The  electric  field  is  reduced  20%  per  decade  of  pressure  between  le-5  and 
1  torr.  This  is  an  exponential  reduction  in  the  field  strength  with  pressure 
and  consistent  with  the  dE/dx  calculation  in  the  text.  This  effect,  coupled 
with  beam  propagation  calculations,  are  well  founded  theoretically  and 
experimental  inferences  are  consistent  with  this  model. 

Neither  direct  observation  of  the  isotropy  of  the  electrons,  nor  the 
change  of  energy  with  pressure  are  available,  since  the  Faraday  cups  were  not 
fielded  on  the  pressure  scans  or  at  varying  azimuthal  angles.  In  future 
work,  placement  of  identical  cups  at  varying  angles  and  maintaining  these 
cups  throughout  the  experiments  is  recommended.  In  addition,  DD2  and  FMM 
data  should  be  passively  electronically  integrated  before  recording  to 
eliminate  digitization  difficulties. 


Chapter  2 


INTRODUCTION 


An  extensive  experimental  session  was  conducted  to  Investigate  the 
behavior  of  Laser  Induced  ElectroMagnetlc  Pulses  (LIEMP)  parametrically 
with  CO2  laser  Intensity  within  a  decade  of  10  GW/cm^.  Parametric 
variations  were  also  made  with  background  air  pressure  from  1  mlcrotorr  to 
a  torr.  Twelve  hundred  (1200)  data  shots  were  taken  with  measurements  of 
electric  fields  on  a  ground  plane  outside  a  glass  vacuum  chamber,  Faraday 
cup  measurements  within  the  vacuum  chamber  30  cm  from  the  Copper  target 
and  30  degrees  off  the  Incident  laser  vector  and  current  measurements  in 
the  target  plane  were  made  with  a  Rogowski  style  monitor  within  a  3  cm 
diameter  of  the  target.  Data  was  recorded  on  oscilloscopes  and 
subsequently  digitized  for  Integration  and  other  analysis.  The  overall 
experimental  arrangement  Is  shown  In  Figure  2.1. 

The  experiments  were  conducted  at  Los  Alamos  National  Laboratory 
using  a  1.3  +  .15  joule  10. I  +  1.  nanosecond  CO2  laser  by  AFWAL/FIESL, 
Technology  Scientific  Services  and  LANL  personnel  after  preliminary 
experiments  performed  at  AFWAL  and  after  joint  planning  with  PRC.  PRC 
provided  and  tuned  the  Faraday  cups.  The  Incident  laser  pulse  was 
asymetrlc  with  a  typical  risetime  of  2  ns  and  an  average  full  width  at 
half  maximum  of  10  ns.  Laser  waveforms  were  monitored  with  a  pyroelectric 
detector  and  energy  was  monitored  with  a  calorimeter.  The  Incident 
waveform  is  shown  In  Figure  2.2.  The  laser  spot  size  was  measured  to  be 
.84  +  .05  mm  FWHM  and  Intensity  was  computed  for  each  laser  shot. 


Intensity  variations  were  created  using  neutral  density  filters  before  the 
100  cm  focal  length  lens. 

PRC  was  retained  to  perform  theoretical  analysis  and  deduce  scaling 
relations  from  the  reduced  data  set.  First  the  theoretical  analysis  was 
done  based  on  the  parameter  space  of  the  experiments  using  existing 
literature  and  original  calculations.  The  reduced  data  was  then  provided 
by  AFWAL  and  PRC  analyzed  the  data  for  empirical  relations  and  compared 
the  theory  to  the  data.  Theory  and  experiment  agree  well  and  form  the 
basis  for  an  understanding  of  the  complex  plasma  processes  and  scaling 
laws  to  predict  subsequent  behavior  of  other  systems. 

The  report  is  organized  as  the  analysis  process  occurred.  The 
theory  is  presented  in  Chapter  3,  the  Data  Analysis  in  Chapter  4,  and  the 
Scaling  in  Chapter  5.  Appendices  are  added  for  details  not  appropriate 
for  the  main  report.  A  detailed  description  of  the  experiments  Is  outside 
the  scope  of  this  report  and  interested  parties  are  referred  to 
AFWAL/FIESL,  Lt.  Adam  Bigelow. 


Chapter  3 


THEORETICAL  ANALYSIS 

3.0  INTRODUCTION  TO  THE  THEORETICAL  ANALYSIS 

Since  the  LIEMP  and  fast  particle  detector  signals  are  caused  by 
fast  electrons  from  the  target,  in  section  3.1  we  evaluate  the  threshold 
intensities  and  gradient  scalelengths  required  for  various  processes  known 
to  generate  hot  electrons  when  an  intense  laser  irradiates  a  target.  We 
conclude  from  these  thresholds  that  the  incident  laser  intensity  in  the 
experiment  is  below  the  threshold  for  production  of  hot  electrons  by 
parametric  instabilities  (the  process  of  interest  in  the  higher-intensity 
laser-fusion  experiments).  Even  allowing  for  the  increase  in  intensity  at 
the  critical  surface  due  to  "swelling",  the  intensity  at  the  critical 
surface  is  probably  below  the  instability  thresholds.  But  in  section  3.2 
we  show  that  the  threshold  and  timescale  for  rippling  the  critical  surface 
is  exceeded,  so  that  resonance  absorption  is  possible.  Resonance 
absorption  Is  known  to  generate  hot  electrons. 

Not  all  of  the  fast  electrons  generated  (by  whatever  process)  near 
n  -  n^  escape  very  far  from  the  target.  Even  when  the  target  is  grounded, 
the  population  of  hot  electrons  with  velocities  away  from  the  target 
creates  a  negative  space-charge  cloud,  or  virtual  cathode,  just  in  front 
of  the  target.  That  space  charge  serves  as  a  barrier  for  all  but  the  most 
energetic  of  the  hot  electrons.  The  self-consistent  physics  of  the  space- 
charge  barrier  is  analyzed  In  section  3.3,  and  it  is  found  that  the 


barrier  thickness  is  related  to  the  hot  electron  Debye  length,  i.e. ,  to 
the  density  and  temperature  of  the  hot  electrons. 

The  current  of  escaping  fast  electrons  is  analyzed  in  section  3.4. 
There  is  no  theoretical  model  that  gives  the  1.75-power  scaling  with  laser 
Intensity  Indicated  by  the  current  loop  monitor,  but  the  analysis  does 
lead  to  electric  fields  scaling  as  observed  in  the  experiment.  (It  may 
be,  as  discussed  in  the  data  analysis  of  section  4,  that  the  loop  monitor 
is  shorting.) 

In  the  experiment,  time-resolved  measurements  were  made  of  the 
arrival  of  some  of  the  escaping  fast  electrons  at  a  Faraday  cup  30  cm  from 
the  target.  Since  the  escaping  electrons  have  to  go  over  a  space-charge 
barrier,  they  are  delayed  (as  compared  with  the  transit  times  calculated 
without  barrier).  In  section  3.5  we  show  that  the  amount  of  this  delay 
allows  estimating  whether  the  target  is  in  fact  electrically  grounded  or 
isolated;  the  delay  when  the  target  is  grounded  is  small  (it  is  calculated 
in  appendix  A),  whereas  electrically  floating  targets  give  large  delays 
(calculated  in  Appendix  B). 

The  fast  electrons  that  escape  beyond  the  barrier  form  a  plume.  In 
vacuum  this  plume  spreads  under  the  repulsive  forces  of  its  own  space 
charge  as  it  expands  outward  from  the  target.  This  is  treated  in  section 
3.6.  The  plume  typically  has  spread  from  a  relatively  narrow  cone  near 
the  barrier  to  a  broad  (tens  of  cm)  front  by  the  time  it  reaches  the 
Faraday  cup.  The  spreading  is  less  when  there  is  appreciable  ionization 
of  a  background  gas,  because  background  secondary  electrons  can  move  to 
reduce  the  net  space  charge,  i.e.,  the  self-repulsion. 

Ionization  of  background  has  a  negligible  effect  for  the  10-^  torr 
"vacuum  case,"  but  when  gas  is  introduced  the  plume  is  more  concentrated 


(less  spread  by  Its  space  charge) •  Is  this  concentration  offset  hy 
colllslonal  slowing  of  the  fast  electrons  theaselves?  At  high  pressures 
of  course  the  fast  electrons  are  slowed  and  scattered  and  do  not  even 
reach  the  detector.  For  Intermediate  pressures  some  estimate  of  the 
collision  processes  is  necessary  and  this  is  provided  In  section  1.7. 

In  section  3.8  we  calculate  the  electric  and  magnetic  fields  of  th 
escaping  electron  plume,  and  the  motion  of  secondary  electrons  under  the 
influence  of  these  fields.  What  limits  the  radial  excursion  of  secondary 
electrons  is  the  experiment  timescale.  Since  they  are  slow,  they  don't 
move  very  far  (a  few  cm).  After  the  pulse  of  escaping  fast  electrons  Is 
over  these  secondary  electrons  move  back  into  the  region  where  the  plume 
was,  since  they  are  attracted  there  by  the  Ions  created  by  the  Ionization 
They  overshoot  and  oscillate  about  the  secondary  Ion  cloud.  This  Is  also 
treated  In  section  3.8.  13ie  space  charge  plume  (fast  electrons  ♦ 
secondary  electrons  -  secondary  ions)  may  be  approximately  thought  of  as 
an  antenna,  loaded  hy  the  t lme -dependent  and  space -dependent  charge 
density.  This  antenna  radiates  an  electromagnetic  field  with  frequency 
Inverse  equal  to  the  timescale  of  the  plume.  This  Is  also  treated  in 
section  3.8  and  Is  distinct  from  the  hlgher-f requeory  elect romagnet l <• 
components  radiated  from  time-varying  space  charge  (etc.)  In  the  suh- 
millimeter-size  space-charge  cloud  of  the  potential  harrier  near  the 
target. 

Later,  In  the  chapter  on  scaling  of  data  (Chapter  S)  we  bring 
together  various  of  the  formulas  derived  In  this  chapter  (3)  to  Interpret 
the  theoretical  scaling  of  fast  electron  current,  current  density  at  the 
Faraday  cup,  electromagnetic  signals,  etc.,  in  conjunction  with  the 
observed  experimental  scaling.  It  should  he  noted  that  the  laser- 


wavelength  scaling  of  the  LTEMP  signals  and  of  the  other  measurable 
quantities  Is  coapletely  contained  In  the  dynamics  of  the  electron  plume 
and  space~charge  layer.  Once  the  hot  electron  temperature  Th  and  the 
escaping  hot  electron  current  density  J  are  calculated  from  the  laser 
Intensity  l  and  wavelength  X  as  In  sections  3. 1-3.5,  the  dependence  of  the 
observables  on  I  and  X  is  coapletely  specified,  except  that  the  weak 
dependence  of  the  function  f  on  X  (Eq.  3.52)  Is  not  known. 


3.1 


THRESHOLDS  FOR  VARIOUS  PROCESSES  KNOWN  TO  PRODUCE  FAST  ELECTRONS  AT 


LASER  TARGETS 


3.1.1  Density-Gradient  Thresholds  for  Instabilities  Generating  Plasma 
Waves 

The  ablation  of  plasma  from  a  laser-irradiated  target  In  vacuum 
produces  a  region  of  nearly  uniform  electron  temperature  with  a  time- 
dependent  density  profile  approximated  by 


exi 


X  +  c  t 


(3.1) 


where  ng  Is  the  density  of  the  solid,  cg  is  the  sound  speed  in  the  plasma, 
and  x  Is  the  distance  from  the  solid  surface.  For  a  CO2  laser 

19  —3 

(A  -  10.6  inn)  the  critical  surface  where  u>  -  u  occurs  at  n  -  10  cm 

P  e 

well  out  In  the  ablation  region  at  all  but  the  earliest  times.  The 
electron  density  gradient  scale  height 


l 


(-1 


3n  -1 
n 

e 


(3.2) 


Is  then,  on  the  spatial  average,  roughly  of  order 


*  -  cgt  ,  (3.3) 

which  Increases  from  £  10  Pm  at  t  -  1  ns  to  £  100  urn  at  t  -  10  ns.  At  the 
laser  Intensities  of  the  experiment, 


I  ~  1010  -  2  *  10llW/cm2  , 


12 


7-10  eV, 


Che  plasma  electron  temperature  is  expected  to  be  of  order  T£  * 
with  a  coexisting  hot-electron  component  of  much  lower  density  and  much 
higher  temperature  Tj,* 

At  least  two  instabilities  are  known  to  occur  at  sufficiently  high 
laser  intensities,  which  can  accelerate  electrons  to  very  high  velocities, 
comparable  to  the  phase  velocities  of  the  associated  waves.  We  examine 
the  threshold  intensities  of  the  ones  with  growth  times  shorter  than  the 
~10  ns  laser  pulse  of  the  experiment.  In  addition,  resonance  absorption 
is  known  to  produce  hot  electrons. 

The  instability  with  lowest  threshold  is  the  two-plasmon  decay,  or 

2up  instability,  in  which  the  laser  electromagnetic  wave  decays  into  two 

electrostatic  plasma  waves  in  the  vicinity  of  o  ■  2u^,  i.e. ,  near  the 

"quarter  critical"  surface,  where  n  -  1/4  n  (Liu  and  Rosenbluth,  1976; 

e  c 

Simon  et  al.,  1983). 

A  second  instability  is  that  of  stimulated  Raman  scattering,  in 
which  the  laser  electromagnetic  wave  gives  rise  to  a  scattered 
electromagnetic  wave  and  an  electrostatic  plasma  wave,  at  a  range  of 
densities  n  <  n£  (Liu,  Rosenbluth,  and  White,  1974). 

In  both  cases,  the  plasma  waves  can  grow  to  sufficient  amplitudes 
to  trap  thermal  ("cold")  electrons  and  accelerate  them  to  velocities  much 
greater  than  thermal  (Manheiraer  and  Klein,  1974). 

The  threshold  for  the  2 instability  in  a  density  gradient  with 
scalehelght  l  is  given  roughly  by 

v  1/2 

ir  '  2  (a)  <3'4) 

Th 


with 


,»  ,*  ,1  >.« 


-2^  -  270  / I(W/c»2 ) (X/ 10 .6  um) 


(3.5) 


and  v_.  ■  A. 2  x  10  (eV).  The  threshold  then  can  be  written  as 

in  c 


I(W/cm2)(X/10.6  Wm)  -  1.6  *  10UT  (eV)/l(um)  . 

c 


(3.6) 


Both  Tc  and  i  Increase  with  time  during  the  laser  pulse,  but  for 


T  ~  8  eV  and  t  ~  85  Um  one  can  see  that  the  threshold  COj  laser  intensity 
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would  be  of  order  1.5  *  10  W/cm  .  We  use  this  example  because  a  soft 


threshold  behavior  for  significant  LIEMP  signal  is  seen  in  the  experiment 


10  2 

at  vacuum  Intensity  of  order  1.5  *  10  W/cm  ,  and  T  is  expected  to  be 


7-10  eV. 


The  threshold  for  the  stimulated  Raman  scattering  (SRS)  instability 


in  a  density  gradient  of  scalehelght  l  is  (Estabrook  et  al. ,  1980) 


*osc  s  ,  /  , c  .2/3 

—  <sr> 


(3.7) 


i«G«  ) 


Ken.)  >  10  (— rJ - y>~3/4(-r<Y  T  um)"1/2 

10l0W/cm2  10,6  Um 


(3.8) 


This  Instability  is  thus  quenched  by  the  short  scaleheight/low  Intensity 


in  the  experiment. 


3.1.2  Hot-Electron  Temperatures  from  Resonance  Absorption 


Even  in  the  absence  of  unstably  generated  plasma  waves,  it  is  well 


known  that  resonance  absorption  near  the  critical  density  can  generate  hot 


•  I’.  t* 


electrons  (Forslund  et  al. ,  1977;  Estabrook  and  Kruer,  1978).  Here  a 
P-polarlsed  wave  Incident  at  a  slight  angle  from  the  local  normal  suffers 
a  resonance  Increase  In  the  longitudinal  E  field  of  the  wave  near  the 


critical  surface,  corresponding  to  Landau -damped  space-charge  waves. 
Quiver  motion  of  electrons  near  the  resonance  leads  to  a  drift  and 
acceleration  of  the  electrons  due  to  ponderomotlve  (wave  pressure)  forces 
(Hora,  1979).  See  Fig.  3.1. 

Several  estimates  have  been  made  of  how  the  resulting  hot-electron 
temperature,  T^,  should  scale  with  intensity  and  wavelength.  The  first 
one  considered  is  (Forslund  et  al. ,  1977) 


\  '  70Toin[(3Tr-T>'^7-’2>‘/3  • 

10  W/cm  10.6  Um 


(3.9) 


The  second  (Estabrook  and  Kruer,  1978)  Is  a  curve-fit  to  simulation  data: 


Th  "  47  Tc°'°4  ((-nr1- 

■  *  ^  1  /"i  *'"'»«  V  1  /■»  r 


(3.10) 


10  W/cm  10.6  pm 


for  IX2  values  between  lO1^  and  10* 7  (W/cm^)(um)2.  In  all  cases,  Th  and 


Tc  are  understood  to  be  In  eV. 


A  third  estimate  (Albritton  and  Langdon,  1980)  Is 


T.  -  9  »  10~4(f IX2T  )0,4 

n  c 


7  111/2 

for  flA  <  10  r  (3.U) 

c 


9  *  10_7(fIX2)2/3 


for  fIX2  >  101  1  t3/2  0.12) 


with  f  a  numerical  factor  of  order  0.3  representing  the  conversion  of 
laser  power  to  hot  electrons,  and  with  1  in  W/cra2  and  X  in  um.  The  high- 
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Fig.  3.1  Resonant  matching  absorption.  The  sinusoidal  line 
represents  the  electric  field  E  of  laser  light  oscillating  at 
frequency  «*.  The  sloping  surface  represents  a  density  gradient  in 
the  plasaa  that  increases  toward  the  target.  Because  the  light  is 
oblique  (at  angle  9),  it  is  reflected  when  the  density  of  the  plasaa 
reaches  nc  cos  o.  Oblique  light  also  has  a  component  of  its  field 
perpendicular  to  the  surface,  and  that  part  of  the  field  can  tunnel 
Inward  to  couple  with  the  longitudinal  electrostatic  field  of  the 
plasaa  wave.  This  coupling  occurs  most  efficiently  at  the  resonant 
matching  point,  that  is,  the  point  where  the  plasaa  density  equals 
n,  and  the  plasaa  frequency  u  equals  u.  The  steeper  the  density 
gradient,  the  smaller  the  gap^between  the  reflection  and  resonant 
matching  points  and  the  more  efficient  the  absorption. 


Intensity  portion  of  this  scaling  law  Is  similar  to  the  estimate  of  Max 
(1982) 


30  [f-‘(  ■„“*  ■)( 

101UW/cn> 


abs 


2i  2/3 


10.6  pm 


-V) 


(3.13) 


with  f  a  flux -limit  factor  for  hot  electrons,  of  order  0.1,  and  lat,s  the 
absorbed  Intensity. 

411  of  these  scalings  predict  ~  200  eV  for  CO2  laser  intensity 
11  2 

of  10  W/cra  and  T^  ~  7  eV.  We  note  that  the  mean-free  path  of  an 
electron  of  energy  e  becomes  longer  than  the  density  gradient  scalehelght 
i  when  £(eV)  >  20^* (urn),  and  for  anticipated  100  pm  gradients  this 
indicates  that  a  soft  threshold  for  hot  electron  emission  should  occur 
around  1  (at  n  )  ~  10  W/cm  for  CO2  lasers.  The  observed  escaping  'tail' 
of  this  hot  distribution  should  be  at  energies  several  times  Th  because  of 
the  space-charge  barrier  that  forms  to  control  the  hot-electron  escape 
rate.  Such  circumstances  appear  to  be  more-or-less  as  observed  in  the 
experiment.  These  scalings  would  also  agree  "reasonably"  with  the  fast 
electron  energies  inferred  from  Faraday-cup-arrlval-time  data  in  the  LANL 
experiment. 

[At  early  times,  one  should  note,  the  cold-electron  temperature  T£ 
depends  on  I  in  a  way  not  made  explicit  in  the  scaling  laws  referred  to 
above,  in  that  the  quiver  energy  of  initially  cold  electrons  gives  a 
temperature 


T  (eV)  >  0.37  (  ..  -A -J( - - - )2  (3.14) 

C  10  W/cm  10.6  Pm 


but  by  the  time  a  5  ns  risetime  laser  pulse  has  delivered  appreciable 


power,  Tc  has  generally  reached  values  an  order  of  magnitude  or  more  in 
excess  of  this  initial  quiver  temperature.] 

Resonance  absorption  requires  density  gradients  at  a  slight  angle 
to  the  direction  of  incidence.  Less-than-opt imal  angles  lead  to  lower 
values  of  Th,  and  probably  to  fewer  hot  electrons.  The  steepening  and 
corrugation  of  electron  density  contours  by  the  laser  probably  is 
responsible  for  the  generation  of  non-normal  incidence  from  a  laser-target 
geometry  in  which  the  laser  is  initially  incident  normal  to  the  target. 
Corrugation  of  density  contours  proceeds  together  with  f ilaroentat ion  of 
the  laser  beam  and  its  local  enhancement  of  peak-intensity  "hot  spots"  in 
the  laser  irradiance  profile  (Estabrook,  1976),  so  that  the  peak  intensity 
at  nc  may  be  larger  than  the  Intensity  in  vacuum.  Filamentary  intensity 
enhancement  and  rippling  of  the  critical  surface  are  discussed  next. 


3.2 


FILAMENTATION :  MECHANISMS,  THRESHOLDS,  AMPLITUDES,  AND  EFFECTS 


Critical -density  surface  corrugations  at  a  copper  target  can  occur 

on  a  time  scale  much  faster  than  ion  motion  because  one  may  think  of  the 

electron  density  as  ng  *  z(Tc,T^)n^;  cold  or  hot  electron  temperatures  Tc 

or  T^  can  often  change  faster  due  to  instabilities  than  can  n^  (Trlpathi, 

Ottinger,  and  Guillory,  1983).  When  the  cold-electron  mean-free  path  near 

the  critical  surface  is  shorter  than  about  1/10,  this  thermal  instability 

mechanism  is  more  important  than  ponderoraotively-driven  f ilamentat ion 

(Kruer,  1985). For  a  7  eV  plasma  at  n  ■  n  *  lO^cm  the  thermal 

e  c 

electron  mean-free  path  is  about  0.1  um. 

The  thresholds  for  thermal  and  ponderoraotive  f ilamentat ion  of  the 
laser  beam  are,  respectively  [Kruer,  1985] 


and 


1010W/cra 


2  >  0,2  ^n  ^  ^10  eV 


T  5  X_  2 

~ )  <!*>  i2 

z 


I010U/cra2 


°-2  TeV  (- 


100  yrair 10. 6  bra^ 


*)( 


0 


(3.15) 


(3.16) 


Both  thresholds  are  exceeded  for  all  intensities  in  the  experiment.  To 
demonstrate  f ilamentat ion  it  suffices  to  estimate  the  e-folding  length  for 
the  weaker  process  and  show  that  it  Is  already  adequate  to  yield  order-of- 


*The  reader  is  cautioned,  however,  that  the  calculation  on  which  this 
conclusion  was  based  ignored  radiative  heat  transport  compared  with 
conduction.  Our  experience  with  simulation  of  comparable  plasmas 
indicates  that  radiative  heat  transport  may  dominate,  and  tends  to  smooth 
out  temperature  perturbations.  Hence,  thermally-driven  f ilamentat ion  and 
ponderomot ive ly  driven  f ilamentation  may  be  of  comparable  importance. 


magnitude  transverse  modulation  of  the  laser  Intensity 


An  estimate  for  the  gain  length  (e-foldlng  length)  of  laser  beam 
ponderomotive  f llamentatlon  near  ng  ■  n£  can  be  made  from  the  review  by 
Palmer  [Palmer,  1972] : 


g  *  (urn)  - 


10.6  wmi 


(3.17) 


with  the  fastest  growing  filaments  there  having  transverse  dimension 


k"*1  (Pm)  ~  5  T^2 


^  10^W/cm2j  ^2^  10.6  um^ 


(3.18) 
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One  can  see  that  for  expected  parameters  (TgV  ~  10  eV,  I  >  10  W/cm  , 
f  ~  100  Urn  density  scaleheight)  the  f i lament at  ion,  on  transverse  scales  of 
order  15  Um  or  less,  can  e-fold  several  times  in  the  course  of  getting  to 
the  critical  density.  If  the  'unperturbed'  vacuum  profile  of  the  laser 
intensity  has  peaks  and  valleys  of  order  10Z,  these  may  totally  filament 
the  laser  beam  and  easily  lead  to  at  least  a  doubling  of  the  peak 
intensity,  quite  independent  of  the  resonance  'swelling'  of  intensity. 

During  the  laser  pulse  (~ 10  s),  copper  ions  moving  transversely  at 

the  sound  speed  can  move  a  distance  of  order 


d(ura)  <  16  T1^2  , 

~  ev 


which  is  consistent  with  plasma  f il amentat ion  on  the  same  scale  as  the 
filaments  of  laser  intensity. 


3.3 


SPACE-CHARGE  LAYER  NEAR  AN  ELECTRON-EMITTING  TARGET  IN  VACUUM 


The  potential  energy  of  an  electron  in  a  space-charge  potential 

2 

♦(r)  la  U(r)  ■  -e4(r).  Since  the  Hamimltonian  1/2  mv  +  U  is  a  constant 
of  motion,  a  two -temperature  Maxwellian  distribution  of  electrons  (nc  » 

-3 

cold  electron  density,  cm  ,  n^  ■  hot  electron  density,  nc  +  “  n)  ^as 

the  form 


-  ( l  mv2  +  U)/T 

a  \  (  ®  ^3/2  2  c 

f(vz’vi)  “  (wr>  nc  e 


.  ■>/■)  -(»-  mv2  +  U)/T 

.  r  m  ^3/2  2  h 

+  'stH  °h  ' 

n 


(3.19) 


2  2  2 

Since  v  “  vz  +  vi»  the  distribution  remains  isotropic  at  all  z  where  the 
Maxwellian  model  is  valid.  At  point  r  where  U  -  U(r),  the  cold  and  hot 
densities  are  expressed  in  terms  of  densities  at  U  ■  0  (taken  at 
r  -  0)  by 

-U(r)/T  -U(r)/T. 

nc(r)  *  nc0e  »  nh<r)  "  nh0  e  *  (3,20) 


We  will  subsequently  derive  a  correction  to  n^  due  to  hot  electron  loss. 
Since  we  are  Interested  mostly  in  values  of  U  »  Tc,  we  will  ignore  the 
cold  electron  density  everywhere  except  very  near  r  ■  0,  where  it  is 
mostly  neutralized  by  ions  anyway. 

Poisson's  equation, 

V2$  ■  +  4*e(n  -  n .  ) 

e  1 


21 


can  be  written  in  terms  of  scaled  variables,  In  spherical  coordinates 


♦  -  e*/T.  -  -U(r)/T  n  /n  -  e 

n  n  hen 

«  •  r/iD  •  *D  ■  <Tb/4">n0'2)‘/2  •  VTh  "  eT 


1_  9_ 
2  3x 


(x 


2  3*. 


3  3il) 

TS  (sine  ijg-) 


* 

e 


(3,21) 


where  we  have  Ignored  the  cold  charge  density  since  it  is  very  well  offset 
by  positive  ions,  whereas  n^  is  not.  For  our  problem,  \)i  <  0  everywhere. 
This  model  is  valid  from  z  -  0  (9  -  */2),  which  we  take  as  a  ground  plane 
(♦  ■  0),  out  to  the  first  minimum  of  4*  (maximum  of  U);  beyond  that  we  will 
assume  only  outward-bound  hot  electrons  and  will  use  a  different  model  for 
n(<>)  as  a  result. 

The  transit  time  of  an  electron  emitted  with  normal  velocity  vq  at 
the  source  z  **  0,  U  »  0,  is 


(3.22) 


where  u(r)  ■  2U(r)/m  •  -  2e$(r,9)/m  evaluated  in  the  direction  (0)  of 
travel;  it  is  assumed  that  for  the  fast  electrons  of  interest  9  does  not 
change  appreciably  during  the  motion. 

IJ(r)  first  increases  from  zero,  and  then  decreases,  i.e.  ,  the 
electron  is  retarded  at  first  by  the  hot  electron  space-charge  cloud 


centered  at  r^  >  0,  then  accelerated  outward  once  it  has  passed  through 


the  cloud. 


Fig.  3.2  Potential  vs  distance  froa  ground  plane,  in  direction  of 

detector. 


Electrons  born  with  vq  greater  than  the  escape  velocity  vf , 


v 


2 

c 


■  aax  u  , 

r 


(3.23) 


reach  a  detector  at  large  r,  while  chose  with  vQ  <  v£  do  not.  The  transit 
time  of  chose  that  do  reach  the  detector  Is  a  function  of  Vq,  obviously. 
The  electron  phase  space  Is  shown  In  Fig.  3.3,  along  with  f(vr)  at  various 


Fig.  3.3.  Electron  phase  space  and  radial  velocity  distribution 


Since  some  of  the  hot  electrons  do  escape,  and  these  are  balanced 
In  the  large"  by  an  Influx  of  cold  electrons  along  the  ground  plane,  the 


potential  at  large  distances  r  »  rq  la  that  of  an  electric  dipole:  the 

negative  charge  -Q  at  z  ■  ♦  r  and  the  positive  Image  charge  *<)  at 

d 

a  »  -r  ,  resulting  In 

q 


!  2Qrn 

i  4  (r  »  r  )  - - r3-  cos  6  .  0.24) 

j  q  r2 

i  If  the  Inductance  or  capacitance  of  the  t arget -grouudpl ane  system 

i 

prevents  or  slows  the  influx  of  cold  electrons,  the  potential  also  can 
have  a  monopole  -Q/r  contribution  (which  goes  over  to  a  form  like  -Q/rq 
i  near  r  ■  0). 

From  Flg.3.3§one  can  see  that  the  fast  electron  density  at  r  <  rq 
la  given  by 

m 

f  f  ( v,r  )d2v.dv  , 

,  ,  1  r 

v  --v  .  (r ) 
r  1 

2.2 

with  Vj(r)  -  wc  ~  u(r)  and  u(r)  =  2U(r)/m,  rather  than  by 


I  f(v,r)d  v|dyr  “  n 


HKrKT 


hO 


This  results  In  a  correction: 


-U(r)/T  /«v  / 2  -  l '  ( r 1 

W  [|  -  J  .etc  IJ-— f - 

1  h 


( e  r  f  c  z  =  — -  f  e  <lw  )  , 

z 


r  <  r 


1  •  ®  •  • 


tfhtrt 


nh0«^(l  - -j  «rfc  (/♦  -  , 


♦  <  ♦  <  0  . 

■  - 


0.25) 


We  note  for  leter  reference  that  erfc(O)  »  1 ,  so  that  it  ♦  ■  ^  Eq.  (3.25) 
gives  a  factor  of  2  reduction  In  n^,  as  can  be  expected  intuitively  froa 
Fig.  3.3.  At  ♦  ■  1/2  erfc  <  0.48,  and  at  +  ■  0  if  >  3  the  erfc  is 

<  .085,  i.e.,  a  correction  of  <  4.3Z  to  the  uncorrected  n^. 

For  r  >  r- 


vf-+v2(r) 


f(v,r)d  vxdvr  , 


v2(r)  =  v*  -  u(r)  . 


Continuity  (along  with  the  quaslstatlc  assumption  3n/3t  *  0)  gives 


r  n<v  > 
r 


n,  <v  >.  r 
I  r  l  q 


(3.26) 


In  spherical  geoaatry,  where  <  >  aeans  an  average  over  f(v,r)  at  r,  i.e.. 


r  *  *  2 

<v  >  -  /  f(v,r)v  dv  d  v.  , 

r  r  r  i 


and  subscript  1  refers  to  r  ■  Tq.  One  can  show  that  for  Maxwellian  hot 
electrons  at  r  -  0,  this  results  In 


—  | A  U ( r ) | /T  r  2 

n  (r  >  r  )  -  n  (r  -  r  )  e  •  (-3-)  , 

h  q  h  q  r 


0.27) 
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where  |AU|  Is  Che  change  In  potential  between  and  r. 

[AU  -  U(r)  -  U(r  )  is  assumed  to  be  negative,  i.e. ,  the  electrons  are 
<1 

accelerated  outward, ]  Thus  in  this  exterior  region  one  has  in  dimension¬ 
less  variables 


3  r  2  1 _  .  3<K  Xq  2*(rq)-* 

+  7m  <8in0  "  T  e 


57  *  sins' 


(3.28) 


since  nh(rq)  -  1/2  nQhe 


*(r  J 


[♦(r  )  <  0],  At  large  r,  of  course,  <l>  *  0 


but  exp(2i|»(r  ))  is  small  so  that  v  ♦  "  0  as  expected. 

q 

To  estimate  the  magnitude  of  ij»(r  )  and  the  location  rq>  we  consider 
the  1-0  analog  of  the  problem,  with  the  approximate  (e^)  form  of  n: 


x  <  x 


(3.29) 


with  boundary  condition  ^(0)  -  0.  The  behavior  of  i|>  at  large  x  will  be 
incorrect  because  of  both  the  1-D  model  and  the  altered  form  of  the  right- 
hand  side.  We  note  that  a  one-parameter  family  of  first  integrals  is 


W  -  C)  , 


+(x) 


0  <  C  <  e 


+(xq)  i  ♦  £  0  . 


(3.30) 


yi|i  \U  7 

e  -  C,  we  have  e  ■  w  +  C  and 


(3.31) 


/T  dw 


-  w4  +  C  , 


Since  w(x*0)  -  » 1  -  C  >  0,  the  solution  is 


V.  AV.V.  aV.V.V.  ,’,V,V 


-  /c  can  (k  -  ax)  for  ex  <  x  , 


(3.32) 


where 


Then  since 


<e  -  can 


2  *  „ 
w  -  e  -  C, 


_1(  J--~-  )  and 


-  c7T 


(3.33) 


we  have 


♦  ■  InC  -  2£n  cos(k  -  ox) 


x  <  x  "  </a  . 

-  q 


(3.34) 


When  x  -  </o  ehe  In  term  is  xero  and  ♦(r  )  -  tnC  <  0.  When  x  -  0, 

q 


Che  condlClon  ♦  ■  0  requires  cos  k  ■  C,  which  la  ldenClcally  true.  The 


K 


The  position  is  then  /2/M  cos  /TT,  where  N  -  C  •  n(x  )/n(0)  <  1. 
For  N  «  1,  cos  /¥  *  x/2  and  r^  *  /2/N  X^,  where  X^  is  the  Debye  length 
of  hot  electrons  at  x  ■  0,  as  defined  previously.  The  potential  barrier 
height  is 


|e*|  -  Thln(l/N) 


at 


r  -  /T7N  Xn  (3.36) 

q  D 

for  the  grounded  target  plane  in  1-D.  To  the  extent  that  r^  is  less  than 
or  comparable  with  the  spot  size,  the  1-D  analysis  should  be  rather 
realistic  for  r  <  1.5  r  . 

q 

Note  that  when  the  potential  barrier  height  is  a  few  times  Th  the 
number  density  of  hot  electrons  at  and  beyond  the  barrier  is  quite  small: 

-|e*|/Th 

N  -  n.  (r  )/n(0)  ~  e  «  1  ,  (3.37) 

h  q 


and  their  mean  energy  at  large  distances,  assuming  4*  +  0,  is  of  order 
<T  in( 1/N)>. 

n 

We  note  in  passing  that  the  dynamics  of  outbound  fast  electrons  is 
nearly  collisionless  near  and  outside  the  critical  surface,  since  their 
mean  free  path  is 


Xmfp(cm)  ~ 


T*/n 


we  (1  +  Z)inA 


[T  (keV)]2  n 

n  C 

[1  +  Z(T)]  n 


1.6 


(3.38) 


for  CO.  laser  n  .  This  Is  easily  seen  Co  be  longer  chan  a  100  urn  density 
gradient  scalehelghC  once  Che  Cyplcal  hoC-elecCron  energy  exceeds  abouC 
160  eV  (for  Z  -  3).  Hoc  elecCrons  Chen  are  decoupled  from  Che  plasma  when 
ouCward  bound,  excepC  for  elecCrosCadc  retardation;  only  Che  inbound  ones 
lose  their  energy  colllsionally  Co  Che  dense  cold  electrons  and  to 
inelastic  processes. 

The  fraction  of  hoc  electrons  near  Che  critical  surface  which 
escape  over  the  potential  barrier  depends  on: 

(a)  their  probability  of  collision  before  escape 

(b)  the  ratio  of  potential  barrier  height  to  T^ 

(c)  oscillations  of  the  potential  barrier 

(d)  the  angular  distribution  of  their  production. 

The  last  of  these  depends  on  the  degree  of  corrugation  of  the 
density  contours  near  the  production  region.  Many  hot  electrons  may 
'reflex'  through  the  strong-field  region,  l.e. ,  be  reflected  by  the 
potential  barrier,  scatter  in  angle  in  the  denser  plasma,  and  be  re- 
accelerated  on  a  second  passage.  About  half  of  the  hot  electrons  produced 
near  nc  probably  enter  the  dense  plasma  and/or  solid  at  n  <  n^  and  lose 
most  of  their  energy  colllsionally  to  the  plasma,  increasing  Tc.  The 
probability  of  momentum  transfer  collision  of  a  fast  electron  with  birth 
energy  traveling  out  from  the  source  in  direction  0  in  a  distance  z  Is 
approximately 


[1  +  Z(T  )]sec0  z  ,  .  T.  2 

PA1»(z’e*0)  " - 7 -  J  dz  4^-  h — rWJ  .  <3-39> 

&ra  1.6  T?  0  nc  E0z  U(Z> 


when  this  is  <  l,  from  Eq.  (3.38).  Here  Gqz  =  e^cos  The  minimum  of 


£A_  -  U(z)  occurs  In  front  of  the  plasma,  where  n(z)  «  n  .  Near  n  , 

UZ  c  c 

2 

where  U/T^  is  small,  the  square  bracket  is  nearly  T^/e^cos  0  and  one  thus 


PAm<8,e0’9) 


1  +  Z(T  )  c  z  .  v 

- sec50  /  dz^- 

1.6  0 


when  <  1  (3.40) 


with  in  keV.  The  free-escape  probability  can  then  be  estimated 
approximately  as 


f>esc<z ,£0’6^  '  1  -PA„(2-Ve)  •  (3-“> 


The  Integral  of  n/nc  is  approximately  the  density  scaleheight. 

Because  of  the  strong  0  dependence,  one  can  estimate  the  0-width  of 


escaping  fast  electrons  of  energy  £q: 


cos  A0(e  )  ~  1.2[£(cm)]0,2[e  (keV)]-0,4  .  (3.42) 


(We  have  taken  Z  ■  3  as  expected  In  the  copper  blow-off  at  7-10  eV. )  For 
£q  ~  T^  no  particles  escape  over  the  barrier  because  U  Z  3  T^.  For 
)  4  T^  ~  1  keV  and  l  ■-  100  pm,  this  expression  Implies  a  plume  angle  of 
order  60  *  [c^(keV)]^*^2  and  an  escaping  fraction 


nf./nh 


i-  exp  (-  U  /T.  ) 
6  max  h 


(3.43) 


for  an  isotropically  produced  Maxwellian.  The  barrier  spreads  the  angular 
width  of  the  electron  distribution  passing  over  It,  but  recontracts  It  on 
the  outward-acceleration  side. 


Fluctuations  In  the  barrier  height  typically  occur  on  timescales 

2*/u  (hot),  which  can  be  written  In  terms  of  the  escaping  current 
P 

approximately  as 

T  (ns)  ~  °*35  an(mm)[e(keV)]1/4tI.(A)]_1/2  ,  (3.44) 

08  C  U  D 

where  Is  the  radius  of  the  space  charge  cloud,  of  the  same  order  as  the 

radius  of  the  laser  spot.  The  magnitude  of  these  fluctuations  cannot 

easily  be  assessed,  but  because  U  /T  is  typically  around  3  or  4 ,  the 

m  h 

ratio  of  current  fluctuations  to  steady  current  Is  around  3  or  4  times  as 
large  as  6U/U  and  could  thus  be  significant.  These  current  fluctuations 
on  timescales  of  order  1/2  ns  can  lead  to  radiated  (non-space-charge) 
E-flelds  with  correspondingly  high  frequencies,  superimposed  on  the  fields 
with  pulse-duration  timescales.  This  will  be  discussed  further  In  section 
3.8. 

The  foregoing  calculation  of  the  space-charge  barrier  Is 
appropriate  In  vacuum,  where  Ionizing  ambient  gas  cannot  provide  charge 
neutralization  by  driving  out  newborn  plasma  electrons.  Except  at  very 
high  pressures  and  hot-electron  current  densities,  Ionization  by  hot 
electrons  Is  probably  negligible,  as  considered  further  In  section  3.7. 

But  photoionization  from  the  heated  target  is  not  necessarily  negligible 

_3 

at  pressures  of  order  10  torr  and  above.  Using  the  approximate  formula 
for  photolonlzat Ion  rate 

(3.43) 

with  the  radiated  blackbody  Intensity 


nnI.(W/cm2)  a  (cm2) 
0  9  ✓  91 

1.6*  Hf1*  T^rr 
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T_  .4  rrt  2 


%(W/o»2)  -  2.5  x  ,o6  t-r-Sjj)  (j^)  «  (3.46) 


in  a  blackbody  spectrum  with  emlssivity  e  and  photon  energies  peaking  at 


~  3  T^,  and  with  photoionization  cross  section  of  the  gas 


-17  2 

o,  ~  2  x  10  cm  , 


(3.47) 


and  gas  density  nn  *  3.5  x  10  P(torr),  we  arrive  at 


T  3 


"i  ~  (1°10»'3  ( 

r  (cm  ) 


-3 

3  x  io  torr 


)  (3.48) 


for  r  »  Tq,  the  Irradiated  target  spot  size. 


At  3  x  10  J  torr,  this  could  lead  to  a  ball  of  ionization  (falling 


off  in  density  at  1/r  )  comparable  to  the  hot  electron  density 


-'Kr)  j 

n  (r)  ~  n  (0)e  q  (r  Ivy 
n  n  q 


after  times  of  order  5  ns,  if  n,  (O)r^  <  3  *  lO1^  cm  *  (for  *  »  4). 

h  q  ~  m 


18  —3  _o 

In  fact,  however,  we  estimate  n,  (0)  ~  10  cm  and  r  ~  10  cm,  which 

h  q 


would  imply  photoionization  density  of  order  IX  of  the  fast  electron  plume 


density  at  small  distances  (10~^  to  10  cm). 


3.4 


CURRENT  DENSITY  OF  ESCAPING  FAST  ELECTRONS 


If  one  assumed  a  current  density  of  fast  electrons  (at  the  target) 

related  to  the  hot-electron  temperature  by  the  Richardson  equation,  J 
2 

would  scale  as  T.  : 

n 

2  10  2 

J(A/cra  )  -  1.6  x  10  U[T.  (eV)]exp[-  e*  /T.  J  .  (3.49) 

n  m  n 

[For  T^  scaling  as  1*/^  (gq.  3.9)  this  would  give  a  frequency  for  the 
virtual-cathode-oscillation  component  of  the  LIEMP  signal  scaling  as  I  '  , 
for  the  vacuum  case.]  But,  in  fact,  the  current  densities  from  electron 
acceleration  at  the  critical  surface  must  be  less  than  the  cold  electron 
source  rate,  rather  than  having  their  saturated  values.  Their  current  is 
given  from  the  acceleration  flux  across  the  E-field  layer  at  the  critical 


(We  use  II  to  denote  the  Incident  laser  intensity,  to  avoid  confusion  with 
current#)  Combining  (3.9),  (3.41),  and  (3.42)  gives 

1  1/8 

J(A/cm2)  <  7.1  x  io7  (10‘*  Wm)2  [— . pjt - 7]  exp [-  e*m/Th] 

10  W/cm 

as  an  upper  bound  on  J  due  to  the  limitation  on  electron  supply. 

But  an  even  more  stringent  limitation  comes  from  the  capacitance  of 
the  system,  i.e. ,  from  the  timescale  for  the  ground  plane  to  resupply 
conduction  electrons  to  the  target  as  it  loses  fast  electrons  over  the 
space -charge  barrier. 

Considering  the  path  of  an  accelerated  electron  as  part  of  the 
circuit  by  which  current  is  returned,  one  has 

(a)  the  accelerating  region  at  the  critical  surface,  treated  as  a 
battery, 

(b)  the  nonlinear  capacitance  of  the  space-charge  barrier,  the 
thickness  of  which  is  NA^  (N  is  a  number  of  order  3,  and 

is  the  hot-electron  Debye  length  in  the  layer), 

(c)  the  inductance  of  the  vacuum  path  of  the  escaping  hot 
electrons  in  the  'fountain',  until  they  return  to  the  ground 
plane.  (The  fact  that  v  “  E  in  this  region  makes  it 
inductive. ) 

(d)  the  resistance  of  the  skin-depth  conduction  layer  in  the 
copper  ground  plane. 


See  Figure  3.5 


The  capacitance  is  of  order 


n_  100  eV  1/2  r.  2 

c(pF)  •  185 '-fer5 

h 


(3.53) 


The  inductance  is  of  order 


L(nH)  ■  (path  length,  cm)  . 


The  skin  depth  in  copper  (conductivity  5.88  *  10'  Mh./m)  is  approximately 


-4  Tr  1/2 

4  (cm)  -  3.3  *  10^  (-■  J---)  ' 

2.5  us 


(3.54) 


Fig.  3.5.  Equivalent  circuit  of  hot  electron  paths  from  a  grounded 
target,  showing  capacitance  of  the  hot-electron  space  charge  layer. 
Inductances  of  the  return  paths  of  escaping  hot  electrons,  and 
resistance  of  the  cold  electron  currents  flowing  in  the  ground 
plane. 


where  is  the  laser  pulse  risetime  (t^/2.5  ns  ~  1).  The  resistance  of 
the  copper  out  to  radius  r  is  then 


R 


1  r  dr' 
a  *  WPJ 
r0 


(rg  is  the  radius  of  the  space-charge  cloud,  of  order  1/2  mm);  i.e. , 

R(n)  -  5.2  x  io~3  (?«A-n.s.)1/2  [jL-r  ln  f-1  °*55> 

Tr  r0 

(where  the  square  bracket  is  unity  for  r  =  1  m. ) 

The  circuit  relation  can  be  simplified  by  approximating  the 
distributed  circuit  as  a  lumped  circuit  with  a  typical  path,  i.e., 
assigning  typical  values  to  L  and  R  (guessing  appropriate  path  dimensions 
for  the  hot  electron  plume).  The  analysis  will  turn  out  to  be  somewhat 
insensitive  to  this,  because  the  capacitance  will  he  found  to  dominate 
(though  not  by  a  large  margin).  The  circuit  equation 

LI  +  RI  +  £  /  I  dt  -  V(t)  , 


when,  Fourier  transformed,  becomes 


(-  io>L  +  R  +  ■1--[a)Cr)  I(w) 


V(o>)  , 


(R  -  ~  +  i  Jr)'1  V 


(3.56) 


-1 


with  T  *  03  ■  t  ,  the  laser  risetime.  From  a  typical  20  nH  inductance 


and  a  2.5  ns  risetime,  we  get 


L/t  ~  80, 


(3.57) 
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:/C  ~  14  n  •  ( 


Tr  ulO18  c,-3  Th  V2 


r/TffeV 


(3.  *>8) 


R  ~  4  x 


2.5  nsi 1/2 


1Q--J  q  . 


(3.59) 


When  t/C  dominates,  we  should  expect  capacitively  limited  current 


I  -  CV/r,  or 
c 


XC(A) 


~  72 


10  cm 


which  scales  as  (n  T  )l/  since  e+  is  a  few  time  T.  .  This  results  In 


current  density  at  the  potential  barrier  of 


J  (A/cm  ) 
c 


V  103  ( 


3  f 2.5  ns> 


10  cm 


rOTTeV' 


(3.61 ) 


for  t/C  >  L/t,  i.e.  ,  for  capacitively  limited  current,  and  for 


e$m/T^  "  4.  The  inductance  L  of  an  average  current  path  Is  not  accurately 


known,  and  larger  currents  than  this  capacitive  limit  can  be  drawn  when 


the  risetime  t  Is  tuned  to  </LC,  l.e.,  when  resonance  occurs.  This  is 
r 


nearly  the  case  (probably  within  a  factor  of  2)  In  the  experiment,  but  the 


resonance  should  be  quite  sharp  because  the  resistance  Is  small. 


Note  that  since  e4  is  a  few  times  T.  ,  If  we  assume  the  scaling 
th  n 


Th  .  uL>V" 


as  in  Eq.  (3.9),  the  hot-electron  current  density  at  the  space-charge 


barrier  in  Eq.  (3.51)  would  scale  as 


iXKUAVSjfllMUf  ILkTVUr*  » 


’ll.1 


.%  AA.V'  .SNjS  -.VXV.'.  A 


One  can  sake  an  approxlaate  dlaenslonal  estimate  of  the  hot -elect  ron 


density,  n^,  near  the  critical  surface  (or  quarter-critical  surface  when 
they  are  generated  by  the  2  wp  instability)  as  follows:  the  energy  flux 
of  hot  electrons, 


nhVh 


-  TV2 
"hTh 


is  some  fraction  f  of  the  laser  input  energy  flux,  l,  .  The  fraction  f 
appears,  in  computer  slnulat ions,  not  to  vary  strongly  with  I ^ ,  although 
it  may  increase  with  1.  Thus  if  one  assumes  the  Th  scaling  of  Eq.  0.9), 
one  has 


-  t^/2f< X  )/X  . 


().h2) 


This  gives  a  capaclt lvelv-l isi ted  current  density  at  the  target  scaling  as 


J 

c 


[  f  ( X  )X 


o.hi) 


until  becomes  as  sisa  1 1  as  the  ringing  time,  where  a  sharp  resonance 

asy  occur.  For  shorter  than  this,  the  current  Is  Inductively  (or 
otherwise)  United  and  scales  with  Tr>  The  ^IX  t  tme  scales  as 


which  Is  wo— what  shorter  at  long  wavelengths  and  Is  almost  Independent  of 
Intensity  (« 

It  la  appropriate  to  aentlon  here,  however,  that  the  experimental ly 

observed  scaling  of  fast  electron  current.  If  the  current  loop  monitors 

are  In  fact  measuring  this,  appears  to  scale  as  1^  to  a  power  near  1.75, 

I  /2 

rather  than  as  1^  as  Implied  by  the  capaclt lve ly-1 lmlted  current  model 

of  this  section.  The  electric  fields  measured,  on  the  other  hand,  do 

2/3+ 

scale  as  approximately  as  they  would  based  on  the  capaclt lve -1 lml t 

theory  presented  here.  (The  "V*  Is  accounted  for  by  the  slight  Increase 
of  cold-electron  temperature  In  Eq.  3.9  with  laser  Intensity;  T  *  l*^.) 
There  Is  no  theoretical  basis  for  total  current  proportional  to 

La 

unless  T  «  (rather  than  T  •  l/^+)  and  J  Is  not  capacltlvely 

limited.  In  the  next  section  and  In  section  4,  we  will  argue  from  fast 
electron  transit  times  and  arrival -t lme  distribution  that  Th  appears  to 
scale  much  more  slowly  than  1^. 
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3.5  PAST  ELECTRON  TRANSIT  TIMES  FROM  ISOLATED  VS  GROUNDED  TARGETS 


t 


3.5.1  Mapping  of  Velocity  Distribution  onto  Arrival-Time  Distribution 

t  “  0, 

Fast  electrons  are  born  at  r  ■  0  with  velocity  distribution  f(v^) 

(for  velocities  In  the  direction  of  the  detector).  They  move  In  a 

potential  well  U(r,t),  slowing  down  as  they  go  toward  the  detector.  When 

U(r,t)  has  Its  t-dependence  much  slower  than  an  electron  transit  time, 

there  is  a  well-defined  escape  velocity  v£  .  Electrons  with  Vq  »  v£ 

arrive  first,  and  those  with  v„  «  v  arrive  last,  after  a  transit  time 

u  e 

t  ■  t  (r,v  ),  where  r  Is  the  distance  to  the  detector.  The  dN  particles 
e  c  r 

between  v^  and  v^  +  dvg  arrive  (in  reverse  order)  between  t(vq  +  dv)  and 
t(Vq).  The  rate  of  arrivals  Is  thus 

_dN  _  dN/dvQ  _  f(vQ> 

-dr  -dr/dvQ  ldT/dvQ] 

since  dt/dvQ  is  always  negative  (at  least  when  U  is  nearly  constant).  The 
transit  time  (for  nearly  constant  U)  Is 


t  ■  f  —  with  v^  “  V?  .1  U(r)  ; 
0  v  0  m 


for  a  monotonies lly  time-varying  U  this  generalizes  to 


T 


with 


I  iH. 

m  3r 


evaluated  at  r,t(r)  for  each  r  of  the  Integration. 

If  particles  are  born  at  r  ■  0  during  a  range  of  values  of  t,  their 
arrival  distributions  are  superimposed: 
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n 


■  MUUWHHUUUI 


dN 

dt 


r  !<vo’c 

| dx /dvQ 


T(v0)) 


dt  . 


If  U(r,t)  is  slowly  time  varying,  dt/dVg  must  be  evaluated  at  t.  If 
U(r,t)  varies  on  the  same  timescale  as  t,  more  convoluted  measures  are 
required. 

3.5.2  Target  Potential  and  Fast  Electron  Transit  Times 

In  Appendix  A  we  show  that  for  the  self-consistent  dipole 
potential,  that  is,  when  the  liberated  hot-electron  space  charge  is 
exactly  balanced  by  a  positive  image  charge,  the  unperturbed  transit  time 
of  electrons 


■*0  m  t/vQ  (3.64) 

is  modified  only  by  the  addition  of  terras  of  order  r^/ vq ,  when 
2  2 

Vq/v£  -  1  =  A  is  unity  or  larger  (i.e.,  for  the  faster,  earlier  arrivals). 
Since  rq  «  r  is  expected,  these  corrections  are  a  few  percent  for  the 
early  arrivals.  (rq  is  the  "charge  radius,"  t.e. ,  the  distance  from 
ground  plane  to  the  most  negative  potential.) 

But  if  the  capacitance  and  inductance  of  the  target  and  ground- 
plane  system  are  such  that  the  image  charge  cannot  completely  match  the 
emitted  electron  space  charge,  then  the  target  is  at  a  negative  potential, 
rather  than  zero  as  in  the  dipole  model.  In  that  case  even  relatively 
early  arrivals  will  be  delayed,  as  we  shall  see  presently.  We  refer  to 
this  as  the  raonopole  case  (in  contrast  with  the  dipole  case). 

In  Appendix  B  we  show  that  for  reasonably  early  arrivals  with 
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A  ~  1,  che  unperturbed  transit  tine  In  the  monopole  case  is  multiplied  by 
the  factor  «''!  +  ( I/A)  : 


T(A  >  1) 


^  / 


1  + 


(3.65) 


2  2 

where  again  Vq/v£  -  1  =  A  and  ve  Is  the  escape  velocity.  But  In  the 
monopole  case  ve  Is  determined  by  Che  net  charge  (l.e.  ,  the  well  depth), 
whereas  In  the  dipole  case  ve  was  determined  by  the  self-consistent 
potential  barrier  height  at  a  distance  of  several  Debye  lengths  from  the 
target. 

When  r  Is  written  as  t  •  (r/vg)F(A),  one  has 


dr 

dv 


0 


I_  +  X_  3F  .  ^1°  .  — 

v0  v0  57  "  vj  "  v( 

t 


2r  3JF 

v  3A 
c 


(3.66) 


where  one  always  has  a  monotonlc  F(A)  >  1  and  (3F/3A)  <  0,  since 


Lim  F  ■  1.  In  terms  of  the  arrival  time  distribution  dN/dt ,  It  Is  t 
A+«> 

(rather  than  v^  or  A  or  v£)  that  Is  the  measurable  Independent  variable, 
and  so  t  ■  (r/Vg)F(A)  must  be  solved  for  Vq(t)  and  Eq.  (3.66)  must  be 
expressed  In  terms  of  t.  For  the  very  first  arrivals  (v^  »  v^ ) ,  which 
are  negligible  In  number  for  a  Maxwellian  velocity  distribution,  t  -  r/v 


0 


so 


and 


dt  /dv 


0 


-r/v 


1  2  «_ 

■  7  mv0/Th 


-rT2/r2 


®  2  -2/t 

7  r  T  /Th 


7  r 


(3.67) 


(3.68) 


giving  an  arrival  time  distribution  which  begins  as 


P(T) 


(3.69) 


dN  r  -u/r^ 

"  3t  "  Te 

T 

with 

®  “  T  r2/Th  *  (3.70) 


This  is  non -expandab ly  snail  a*  t  ♦  0  but  becones  significant  at 


T 


(3.71) 


i.e. ,  at  Vq  «  / ft^/m,  if  the  corrections  (?)  to  free-streaming  have  not 
already  become  important  by  the  time  the  arriving  "v^'s"  have  dropped  to 
this  value.  In  the  case  of  the  self-consistent  dipole  potential, 
electrons  with  this  nearly  thermal  speed  are  all  reflected  and  never 
appear  at  the  detector,  so  even  before  time  Tj  the  arrival  time 
distribution  P(r)  begins  to  fall  below  the  unperturbed  values  of  Eq. 
(3.69).  In  the  case  of  the  monopole  potential  due  to  net  charge  Imbalance 
Q  (which  depends  on  circuit  conditions  at  the  ground  plane) ,  thermal-speed 
particles  may  be  transmitted  to  the  detector  if  Q  is  small  enough,  and  if 
no  dipole  potential  barrier  reflects  them.  Generally,  the  highest  barrier 
determines  ve,  a  \d  only  Vq  >  ve  are  detected  (which  implies  that  at  least 
by  time  t  •  r/v£  the  arriving  pulse  has  significant  dispersion  so  that  it 
has  values  below  the  curve  of  Eq.  (3.69)  in  Figure  3.6). 


Fig.  3.6.  Arrival-time  distribution  for  a  Maxwellian  burst 
at  t  •  0,  without  retarding  potentials. 
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The  technique  just  described  for  inverting  the  arrival-time 
distribution  to  Infer  the  energy  distribution  of  the  fast -electron  source 
is  useful  only  when  the  source  duration  is  short  compared  with  the 
arrival -time  distribution.  Wien  both  the  source  strength  and  the  source 
distribution  of  energies  vary  with  time  on  the  same  timescale  as  the 
arrival -time  distribution,  not  enough  information  is  contained  in  the 
latter  to  infer  both  the  source  strength  and  the  temperature;  the  mapping 
is  a  many-to-one  mapping. 

When  the  source  energy  distribution  can  be  approximated  as  fixed, 
however,  one  can  again  in  principle  invert  the  arrival  time  distribution 
to  infer  source  strength  provided  the  delay  introduced  by  the  potential 
barrier  is  small;  each  arrival  moment  t  samples  velocity  v  of  the 
distribution  weighted  by  the  source  strength  at  t  -  L/V. 

If  the  source  flux  is 

F( v ,t)  -  S(t')  /  f , ( v)v (dv |  (3.72) 

with  f,(v)  the  velocity  distribution  integrated  over  v^  ,  then  the  arrival¬ 
time  distribution  is  proportional  to 

« 

P(t)  -  f  S(t  -  L/V)f  (v)v  dv  .  (3.73) 

0 

Dropping  the  subscript  <  and  transforming  to  the  variable  t  ■  L/v  this 
becomes 

t 

P(t)  -  f  S( t  — r )G( t  )dx  ,  (3.74) 

0 
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where 


G(t)  -  [f(v))|  L2/t3  .  (3.75) 

I v»L/t 

Because  this  Integral  equation  is  a  convolution,  it  can  be  formally 
solved  by  Laplace  transforms.  Denoting  the  transform  of  G(x)  by  G(s), 

A  A 

that  of  S(t)  by  S(s),  and  that  of  P(t)  by  P(s),  we  have 

P(s)  -  S(s)G(s)  (3.76) 


so 


S(t )  -  ^-1(P(s)/G(s)]  ,  (3.77) 

where  £  *  denotes  the  inverse  Laplace  transform  operator. 

Shown  above  the  t  axis  in  Figure  3.7  are  the  orbits  x(t,t')  of 
electrons  leaving  the  target  at  time  t'  and  detected  at  x  -  L  at  times  t 
and  t  +  dt.  Shown  below  the  t  axis  are  the  velocities  which  lead  to 
detection  between  t  and  t  +  dt,  for  electrons  born  at  t' ,  over  a  range  of 
t"  values  from  0  to  t  or  t  +  dt.  If  one  imagines  the  source  velocity- 
times-density  distribution  as  a  product  of  hot  electron  density  nh(t')  at 
the  source  times  a  velocity  distribution  (with  unit  area)  f(t',v),  the 
domain  of  this  function  is  the  lower  half-plane  of  Figure  3.7.  The 
density  of  the  orthogonal  cross-hatching  in  the  figure  represents  the 
height  of  n^(t" )f (t' ,v)  on  the  strip  of  (t,v)  space,  for  a  pulsed  source 
with  a  typical  f  monotone-decreasing  in  v. 


However,  In  modeling  the  integrand  of  Eq,  (3.72)  by  Maxwellian  or 

power-law  velocity  'tails'  with  t  ■  L/v,  one  must  take  heed  that  there  is 

only  a  finite  region  in  (t,t')  space  on  which  the  integrand  is  nonzero. 

The  power-law  energy  distribution,  which  gives  simpler  Laplace  transforms 

than  a  Maxwellian  model,  has  both  a  minimum  and  a  maximum  velocity  for  its 

applicability;  and  fast  electrons  cease  to  originate  from  the  target  very 

soon  after  the  laser  source  is  turned  off.  This  delimits  the  region  of 

nonzero  integrand  as  shown  in  Figure  3.8,  and  the  (infinite-interval  0-“) 

Laplace  transform  is  not  that  of  a  power,  but  of  a  function  which  is  a 

power-law  inside  the  parallelipiped  of  Figure  3.8  and  zero  outside  it. 

This  complicates  solving  the  convolution  integral  equation  of  the  first 

kind  (Eq.  3.74)  by  the  Laplace-transforra  procedure. 

For  t  >  L/v  (as  in  the  region  of  nonzero  P(t)  in  the 
max 

experiment),  the  integral  representing  the  arrival-time  distribution  is 

P(t)  «  /  vn.  (t-r)f(v  -  L/r)dT/T^  ,  (3.78) 

Lt  h 

with  upper  limit  equal  to  the  minimum  of  t  and  L/vmj,n,  and  lower  limit 
equal  to  the  maximum  of  t  -  x  and  L/vmflX.  (See  Figure  3.8.) 

The  formal  Laplace  inversion  of  this  to  give  S(t)  is 

S(t)  -  X,~l  \  —  P(t.)— T]  ,  (3.79) 

«£f(L/T)T 

when  it  exists. 

For  a  "constant"  source  strength  during  the  pulse  duration  t  ,  with 
an  energy  distribution  f(e)  *  e  *  between  emin  and  emax*  one  can  easily 
show  that  the  arrival  time  distribution  P(t)  has  the  form 


t 


Fig.  3.8.  Integration  region  in  t,T  space.  Integration  to  give 
F(t)  covers  the  shaded  are  i,  where  the  solid  vertical  line 
represents  the  instantaneous  value  of  t. 


P(t)  -  Min(t,L/vmln)  -  Max  (t  -  tp,  L/v^) 
shown  in  Figure  3.9. 

In  a  typical  data  example  (Figure  4.4d),  we  have  (L  -  30  cm) 

L/v  ,  »  30  ns,  L/v  ■  20  ns,  t  ■  10  ns.  From  L/v  ,  we  get  a  minimum 
min  max  p  min 

energy  escaping  over  the  potential  barrier  e  ■  280  eV,  and  a  maximum 

min 

energy  e  *  640  eV  corresponding  to  the  first  arrivals.  If  we  assume  a 

(BAX 

power  law  f  «  e  *  from  e  ,  to  e  and  correlate  this  with  the 
~  min  max 

temperature  of  a  Maxwellian  distribution  most  closely  resembling  this 
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435  eV,  baaed  on  Che 


power  law,  l.e. ,  f  *  exp  (-e/T^),  we  get  * 
assumption  of  a  source  strength  uniform  In  time  over  the  pulse  length 
with  approximately  constant  Th. 


Fig. 


-2 

3.9.  P(t)  for  "constant"  source  over  t  and  f  *  v  from 

□ 


3.6 


RADIAL  SPREADING  OF  THE  FAST-ELECTRON  PLUME 


The  radial  expansion  of  the  fast -electron  plume  under  the  Influence 
of  its  self-field  depends  on  the  degree  of  charge  and  current  neutraliza¬ 
tion  by  background  ions  and  electrons,  and  can  be  approximately  described 
by  an  envelope  equation  (Lee  and  Cooper,  1976;  Lawson,  1975)  when  the 
expansion  is  not  so  extreme  as  to  make  >  v^. 

At  very  low  background  pressure,  e.g.,  10-^  torr  (the  "vacuum"  case 
in  the  experiment)  the  space-charge  and  current  of  the  fast  electrons  are 
not  neutralized,  and  electrostatic  repulsion  forces  dominate  over  magnetic 

pinching  forces  by  a  factor  of  c/v  .  For  that  case  the  magnetic  field 

f  e 

can  be  ignored  in  the  beam  envelope  equation.  More  generally,  the  net 
expansion  depends  on  the  parameter 


tc 


2  2 
"W 


mv. 


0  -  f 


2/  h 

-  vb/c  ) 


(3.80) 


where  rg  Is  the  radius  of  the  beam  of  fast  electrons  at  the  space-charge 
virtual  cathode  source,  ng  is  its  density  there,  fe  is  the  fractional 
neutralization  of  space  charge  by  ions,  and  v^  is  the  speed  of  the  fast 
electrons  outward  from  the  source. 

In  terras  of  the  parameter  x,  the  beam  envelope  equation,  neglecting 
any  rotation  and  any  slowing  down  of  the  beam,  is  approximately 


d^r 

dz2 


X 

T 


(3.81) 


For  initial  conditions  dr/dz  =  0  at  z  =  0,  this  has  the  implicit 


_  /tn(r/rQ) 

4  /  exp(t2)dt  .  (3.82) 

0  0 

Numerically,  when  f  «  1  we  find 

✓77?  -  d  lie  V  i  .  (3.83) 

Since  I  ~  3A  and  e  ~  1  keV  typically  In  the  experiment,  / < / 2  Is  typically 
of  order  0.3.  We  show  tabular  values  of  the  Integral  In  Eq.  (3.82)  as  a 
function  of  t/tq  in  the  following: 

r/rQ:  10  50  100  300  103 

/:  4.2  15.35  27.4  70.9  209 


Thus  for  Tq  ■  1  mm,  one  has  a  predicted  plume  radius  of  order  30  cm  at 

z  ■  24  era  when  there  is  negligible  space  charge  neutralization.  Although 

this  large  spreading  violates  the  paraxial  assumption  of  the  envelope 

equation,  we  take  It  as  Indicative  of  the  expected  spread.  In  the  vacuum 

and  short-pulse  limits.  Only  after  Ionization  of  background  neutrals  can 

build  up  f  “1  does  the  remainder  of  the  fast-electron  stream  remain 
e 

focused. 
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v’V 


1.7 


FAST  ELECTRON  SLOWING-DOWN  AND  IMPACT  IONIZATION  IN  LOW-PRESSURE 


GAS 


The  range  of  Che  fast  electrons  In  air  or  N2  Is  roughly 

r  *  7  cm/P(torr)  for  e  -  1  keV 

and 

r  -  42  cm/P(torr)  for  c  -  5  keV  , 

-2 

so  at  p  £  10  Corr  where  Che  observed  EM  fields  are  strongest,  the  fast 
electron  range  Is  longer  than  the  distance  to  the  detectors.  At  the 
higher  pressure  P  ~  1  Torr  (where  the  LIEMP  effect  Is  weak)  the  fast 
electrons  are  stopped  before  reaching  the  Faraday  cup.  This  clearly 
affects  the  high-pressure  scaling  of  the  LIEMP  fields  from  the  fast- 
electron  emission  plume,  but  does  not,  per  se,  alter  the  fields  due  to  the 
localized  hot -electron  space  charge  layer  near  the  target. 

The  Ionization  by  the  fast  electrons  can  be  represented  by 

dn 

IT  ’  "0  !  V0v(t)ffe 

*  4.7  *  10^P(torr)n  (cm  at  1  keV 

f  e 


or 


dni  -1 

(cm  /ns) 


~  1 .54  * 


10*^P(torr)J(A/cm^) 


at  1  keV 


where  J  Is  the  fast  electron  current  density. 

The  tons  so  produced  can  move  only  on  a  slow  (ion  motion) 
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timescale,  while  Che  secondary  electrons  may  be  expelled  to  fairly  large 
(cm  or  meter)  distances  on  Ion  timescales  by  the  (predominantly)  space 
charge  field  of  the  fast  electron  beam.  (At  large  currents,  the  magnetic 
field  of  the  beam  tends  to  limit  the  radial  excursion  of  the 
secondaries.)  The  plume  space  charge  Is  thus  neutralized  on  a  timescale 

TgN(n8)  -  0.2/P(torr) 

when  this  Is  less  than  the  beam  pulse  duration  and  when  the  secondary 
electrons  are  expelled.  For  the  low  pressures  at  which  the  peak  LIEMP 
fields  were  seen,  one  has  >  t ,  i.e.  ,  space  charge  neutralization  Is 
Incomplete.  Once  the  pressure  is  large  enough  for  the  neutralization  to 
take  place  In  ~1  ns,  the  space  charge  electric  field  pulse  Is 
significantly  shortened  and  corresponding  E-f leld-drlven  effects  should 
become  much  less  significant. 


3.8 


PLUME  ELECTROMAGNETIC  FIELDS  AND  THE  MOTION  OF  SECONDARY  ELECTRONS 


IN  THE  COLLISIONLESS  LIMIT 


3.8,1  Motion  of  Secondary  Electrons 

A  staple  picture  of  plume  space  charge  neutralization  in  the  low- 
pressure  regime  assumes  rapid  radial  expulsion  of  the  secondary  electrons. 
As  the  beam  space  charge  becomes  compensated  by  the  positive  ions,  we 
might  ask  whether  the  magnetic  field  of  the  beam  will  be  strong  enough  to 
confine  the  secondary  electrons  inside  some  beam  halo  well  before  complete 
electrical  neutralization  is  obtained,  since  Er  can  be  less  than  B^c.  We 
estimate  the  effects  of  the  magnetic  field  in  this  section;  for  beam 
parameters  of  Interest  we  conclude  that  "free"  escape  of  the  secondaries 
from  the  beam  region  is  a  reasonable  assumption  through  much  of  the  charge 
neutralization  phase  at  pressures  below  the  transition  pressure  estimated 
subsequently. 

At  pressures  well  below  the  transition  pressure,  the  secondary 
electrons  might  accumulate  enough  density  near  the  beam  to  be  a  potential 
worry  as  far  as  two-stream  type  instabilities  are  concerned.  But  typical 
orbits  place  the  average  position  of  these  electrons  far  from  the  beam 
region,  so  they  cannot  contribute  significantly  to  the  cascade  ionization 
inside  the  beam  of  low  pressure. 

We  consider  the  fields  and  secondary  electron  motion  that  would 
result  from  a  Bennett  beam  profile  and  a  constant  beam  radius.  Since  the 
Impact  ionization  by  fast  electrons  gives 


3n  , 


>t  “  VinB  • 


(3.84) 
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(with  rig  the  "bean"  density  of  fast  electrons),  we  have  for  the  fractional 
charge  neutralization 


*8 


V  * 


for  t  <  v 


-1 


(3.85) 


with  vJ  the  primary  Ionization  rate  (v.  -  lOr  In  air). 
1  1  torr 


Photolonlzatlon  from  the  target  emission  has  been  neglected.  If  the  fast- 
electron  beam  current  is  Ijj  and  a  fraction  ffi  of  its  charge  is  neutralized 
by  background  ions,  then  the  fields  are 


-  (1  - 


fe> 


^oV 


2v  v 


(3.86) 


2*  (a  +  r  )  zB 


MoV 


- 2 - T~ 

2r(a-£  +  r  ) 


(3.87) 


We  consider  a  secondary  electron  born  at  rest  at  a  radius  r  *  r^.  The 
motion  is  treated  in  the  collisionless  limit,  the  fields  are  assumed 
constant  in  time  during  the  motion,  and  we  neglect  the  Ez  fields  due  to 
the  t ime -dependence  of  the  neutralization  fraction  and  the  attraction  of 
the  ground  plane.  All  these  neglected  factors  will  generally  enhance  the 
escape  of  the  secondaries  at  low  pressure,  although  collisions  will 
ultimately  retard  the  escape  rate  at  sufficiently  high  pressures  where 
mean-free  paths  become  smaller  than  the  beam  radius. 


The  equations  of  motion  with  v^  =  0  are 


dv 


It- 


~evrB0  * 


(3.88) 


dv. 


17“ 


s(Er  “  VzV  ’ 


(3.89) 
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a 


de 

It 


(3.90) 


«  ev  E  , 
r  r 


with 


1  2  2 
e  -  7  ■<%  +  v‘) 


Using  vr  -  dr/dt,  from  Eqs.  (3.89)  and  (3.90)  we  have 


e/ijc 


I„  2  .  2 

(1  -  f  )  -£__I  tn  r  +  “ 


6  VzB  *0 


T"  I  • 

ri  +  a 


(3.91) 


with  the  initial  condition  v  ■  0  at  r  ■  r^,  and  the  definition 


4*  me 


-  17, 000  amp  s  . 


(3.92) 


From  Eq.  (3.88)  we  have  for  constant  Ig 


[B  f  r2  ♦  a2 
In  -y— -j  . 

0  r,  +  a 


(3.93) 


From  Eq.  (3.89)  we  obtain 


to  - 


f e)<c/v^g) 


.  2  2  2  .  2 
1/T/Tv.  r  +  a  i  .  r  a 

7  (1b/I0)  ln  -T"-V  tn  ~ 


T — r 

ri  +  8 


7—T 

ri  +  * 


(3.94) 


The  solution  for  the  orbit  r(t)  from  Eq.  (3.94)  Is  complicated,  but  we  are 
mainly  Interested  ln  the  outer  limit  of  the  orbit,  r^.  For  constant  IR 
this  is  where  vf  *  0,  which  requires 


*>6 


A 

,s 


2  A  2 
rM  +  ‘ 

2  2 

ri  +  * 


2  V1  ~  fe> 
_IB(v*b/c)  * 


(3.95) 


The  smallest  value  of  rM  la  obtained  with  particles  starting  near  the  axis 
(r^  «  a);  their  turning  radius  la  explicitly  given  by 


2I„(1  -  f  > 


m  exp 
a 


b'VeB 


f-1  -  1  • 


(3.96) 


As  f  ♦  1 ,  rM  ♦  r^  •  0.  But  we  see  that  r*  is  larger  than  a  for 
almost  all  fe,  especially  for  currents  much  less  than  17  kA. 

However,  the  actual  maximum  radial  displacement  of  secondary 
electrons  during  the  time  the  beam  of  primary  fast  electrons  Is  flowing. 
Is  r'  ~  (vr/c)tB,  l.e. , 


rM  ~  <ctb)(IB/15  ’tA)(c/v*B)tn(1  +  r2/a2)  .  (3.97) 


For  3A  of  1  keV  primary  electrons  (v  /c  ~  0.06)  lasting  50  ns,  this  gives 

Z  D 

a  radial  traverse  of  order  5  cm,  which  then  subsequently  reverses 
(secondary  electrons  collapse  back  toward  the  axis)  and  oscillates,  after 
the  primary  fast  electron  current  has  passed. 

The  fast  electron  plume  may  be  thought  of  as  having  set  up  a 
nonlinear  space-charge  perturbation  In  the  plume-created  plasma  of 
secondary  electrons,  with  a  frequency  of  order 


f  -  3  *  108  [P(torr)T  (ns)l  (A)/e  (keV) ] '^2/a(cm)  ,  (3.98) 

sp  D  D  a 


which  Is  usually  comparable  with  the  laser  pulse  length  In  the  experiment 


■ 


and  a  ~  3  cm 


,  _  ”2  torr 

for  P  ~  10 


The  apace  charge  E  field  due  to  the  fast-electron  beam  is 


E  (V/cm)  ■  950  I.( A)/ [a(cm)/e  (lteV) )  at  r  -  a  (peak  E)  (3.99) 

r  b 


for  bare  unneutralized  current  Ijj»  with  radius  a.  In  a  raore-or-less 


cylindrical  beam. 


If  the  space  charge  is  neutralized  by  secondary  electron  expulsion 


at  a  distance  d  ~  a  behind  the  beam  front,  the  field  strength  r  at  r  *  a 


Is  reduced  to  2/3  of  the  above  value.  (The  largest  value  of  Er  occurs  at 


r  -  a. ) 


The  Inductive  Ez  field  due  to  B0  Is  approximately 


E  (V/cm) 

z 


2  *  10  ^  I,  (A.)  — — —  at  r  *  a  , 

D  T 

br 


(3.100) 


where  x  Is  the  beam  net-current  rlsetlme. 
br 


From  these  formulas  one  can  see  that  for  a  10  ns  beam  pulse 


rlsetlme  without  space-charge  or  current  neutralization,  the  space  charge 


electrostatic  field  dominates  over  the  inductive  E  field  for  a  values 


expected  in  the  experiment. 


3.8.2  Fields  of  Time -Dependent  Space-Charge  Pltroe 


When  the  Fourier-analyzed  current  In  a  center-driven  linear 


traveling  wave  antenna  of  half-length  L  is 


Ig(<i))exp(l(jJ  |  z  |  / c )  , 


z  <  L 


(3.101) 


'  VLY.V.V\N '  ■.\\\n  v-v*  V-7- .  ■  .•  vV’ 


and  Che  center  Input  current  is 


IQ(t)  -  (2w)  1  /  Iq(cd )exp( -iut )dtu  ,  (3.102) 

and  when  L/c  Is  small  compared  to  the  pulse  duration  of  the  signal,  the 
radiated  field  is  approximately 

E(r,t)  -  (po/eo)1/2(sin0/2irr)(L/c)  ~  I(t  -  r/c)  (3.103) 

when  there  are  no  reflected  waves  in  the  antenna  Itself,  and 

2 

E(r,t)  *  (un/en)1/2(sin6/2irr)(L/c)2  —■  I(t  -  r/c)  (3.104) 

3t 

when  there  is  a  reflected  wave  of  intensity  equal  to  the  outgoing  wave  on 
the  antenna,  l.e.  ,  when  the  Fourier-analyzed  current  is  a  standing  wave, 

I()(“HexP(iwM /c)  “  exp[iw(2L-|z|  / c ] }  .  (3.105) 

(In  both  cases  $  is  polarized  transverve  to  t  but  in  the  rz  plane.) 
(Sengupta  and  Tai,  1976). 

Between  these  limiting  cases  the  radiated  field  is  more 
complicated.  The  radiated  waveform  is  distorted  from  the  input  driving 
voltage  by  both  the  radiation  mechanism  itself  and  the  reflections  from 
the  ends  of  the  antenna.  One  can  see,  however,  that  transients  in  the 
antenna  current  may,  up  to  a  point,  give  large  oscillatory  signals. 
Resistive  loading  of  the  antenna,  l.e.,  suppression  of  reflected  waves, 
causes  the  radiation  field  of  a  transient  to  decrease  at  late  times. 
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On  the  other  hand,  when  the  current  signal  In  the  antenna 
propagates  much  more  slowly  than  c,  one  can  for  comparison  calculate  the 


quasl-statlc  fields.  In  the  experiment  there  are  both  electrostatic, 
partly  radial  (at  the  receiver)  fields  as  well  as  transverse  induced  fields 
at  the  E  probe  (1  m  from  the  driving  point);  the  detector  Is  polarized  to 
pick  up  the  transverse  field.  Including  the  electrostatic  part  parallel  to 
the  antenna. 

The  electrostatic  part  of  the  E  field  arises  from  the  dlpole-like 
space  charge  distribution  and  may  be  estimated  in  the  quasistatic  limit 
from 


E 

z 


ene(r',z')r'dr'dz' 
y (r-r')2  +  (z-z')2 


(3.106) 


using  a  time-dependent  model  for  nfi(r',z')  representative  of  the  expected 
space-charge  plume.  Including  the  halo  of  expelled  background  electrons 
and  their  remaining  Ions  where  appreciable,  and  Including  the  advance  of 
energetic  ions  at  late  times. 

For  a  uniform  line  charge  density  -p  on  the  axis  from  z  *  0  to 
z  ■  L,  and  Its  Image  charge  +p^  at  z  <  0,  the  electrostatic  Ez  Is 


Ez  "  PL{‘ 


-  2 


J  2  2 

y  r  +  z  -  z 


(7 7 


(L-z)2  +  (L-z)  / r2  +  (L+z)2  -  (L+z) 


(3.107) 


In  the  static  limit  where  p  and  L  are  small.  This  fietd  becomes  dipole- 

la 

like  for  r  »  L.  Modifications  to  this  quasistatic  field  arising  from 
finite  radial  extent  of  the  plume  can  be  calculated  from  an  assumed  plume 


60 


shape  (which  may  also  change  slowly  with  time),  but  the  complication  makes 

the  result  less  useful.  It  is  sufficient  to  point  out  from  this  formula 

that  once  the  plume  length  L  Is  several  centimeters,  the  quaslstatlc  space 

charge  field  at  the  D  detector,  1  m  from  the  source  and  near  the  ground 

plane,  can  give  a  reading  even  though  It  Is  not  strictly  a  radiative 

field.  This  portion  of  the  E  field  follows  changes  In  and  l.  in  the  DC 

manner,  l.e. ,  the  behavior  of  Ez  (in  the  absence  of  fast  transients  in 

Pl,L)  mimics  that  of  changes,  pulsations,  etc.,  in  the  net  charge  at  and 

2  3  2  2 

near  the  plume.  It  is  proportional  to  pT L  /r  for  r  »  L  . 

Nondispersive  modeling  of  the  signal  propagation  at  the  fast 

electron  plume  and  its  secondary  electron  halo  would  be  justified  only  If 

the  excitation  frequency  spectrum  were  to  lie  well  above  the  natural 

resonance  frequencies,  e.g. ,  the  range  of  u  space  charge  oscillations 

P 

that  convert,  e.g.,  at  density  gradients,  to  EM  waves  radiated  outward 

into  regions  where  id  <  id.  Neglecting  collisions,  the  local  dielectric 

P 

function  for  frequency  id  transverse  waves  Is 

2 

id 

eT(ui)  -  1  -  -£■  ,  (3.108) 

ID 


2  2 

where  id  -  4ire  n  (r)/m.  This  just  follows  from  the  colllslonless  response 
p  e 

of  the  background  plasma  electrons 


(3.109) 


(where  •  indicates  3/dt  and  the  v*Vv  term  has  been  neglected)  once  the 

Fourier  transform  of  id^E  has  been  approximated  by  -iiDto^'g. 

P  P 

Somewhat  more  generally,  one  has  before  Fourier  transforming 


(3.110) 


VxVxe  +  -y-  (E  +  u>pE)  -  -  jfe  , 

c  c 

2 

where  may  be  a  function  of  r,  z,  and  t,  and  where  is  the  fast- 

electron  current  density.  Because  of  the  3 z3  terras  the  Er  and  Ez  fields 

are  coupled  in  this  vector-pair  of  partial  differential  equations. 

By  taking  the  divergence  of  this  equation  and  using  the  continuity 

equation  for  n^,  one  can  of  course  show  that  the  net  charge  density 

oscillates  with  frequency  w  : 

P 

(.<*1  ♦  3?)(n  )  -  (v.  /c)3?n.  .  (3.111) 

p  t  ne  t  b  t  b 

After  passage  of  the  fast  electrons,  this  represents  pulsation  of  the 
remaining  column  of  secondary  electrons,  although  one  should  note  that 
fast  ions  follow  the  fast  electron  pulse.  To  the  extent  that  Ez  is 
quaslstatlc  and  mimics  these  oscillations,  its  oscillation  frequency- 
squared  should  scale  linearly  with  pressure,  because  until  collisions 
actually  Impede  the  fast  electrons  (at  P  >  0.1  torr) ,  the  number  density 
of  secondaries  is  linear  in  the  gas  density,  for  a  given  fast-electron 
pulse  profile. 

A  fast  electron  plume  which  has  duration  u»  * ,  with  u>  determined 

P  P 

approximately  from  the  peak  secondary  electron  density  near  the  end  of  the 
pulse,  should  set  up  the  largest  oscillations;  thus  for  any  pressure  there 
is  a  pulse  length  which  maximizes  the  LIEMP  (except  for  high-frequency 
components  due  to  transients).  Conversely,  if  the  pulse  length  is  fixed, 
there  should  be  a  'resonant'  pressure,  once  the  fields  due  to  pulsating 
secondaries  exceed  the  vacuum  field  of  the  fast  electron  pulse  itself. 


At  very  low  pressure  then,  when  secondaries  can  be  neglected,  the 

Ez  fields  should  scale  in  the  same  way  as  the  current  of  fast  electrons 

measured  by  magnetic  loops.  At  higher  pressures,  Ez  should  display  an 

1/2 

oscillating  'tail"  and  its  amplitude  may  scale  as  /t^  (since 
2 

u  *  n  above. ) 
p  net 

3.8.3  Radiation  Field  from  Oscillation  of  the  Space-Charge  Cloud  near  the 
Target 

The  "virtual  cathode"  formed  by  the  fast  electron  reflection  layer 

near  the  target  surface  is  not  a  completely  static  feature,  but  oscillate* 

about  its  equilibrium  configuration  (Birdsall  and  Bridges,  1961),  as  well 

as  changing  in  time  as  the  fast-electron  production  rate  changes.  In 

addition  to  radiation  produced  by  the  deceleration  of  electrons  as  they 

are  returned  to  the  target  by  the  potential  barrier  (Barkhausen  and  Kurz, 

1920),  radiation  is  also  produced  by  the  oscillations  of  the  position  and 

shape  of  the  electron  cloud  (Kwan,  1984).  This  radiation  has  typical 

frequencies  on  the  order  of  (5/2)  w  (Walsh  and  Sullivan,  1985)  (w  being 

P  P 

the  plasma  frequency  corresponding  to  the  fast-electron  current) 

f  (GHz)  -  0.32  [J(A/cm2) 1 1/2  .  (3.112) 

osc 

The  radiation  frequencies  of  the  individual  reflected  electrons  are 

also  of  this  order,  because  the  potential  barrier  scale  size  over  which 

the  reflection  takes  place  is  of  order  3  the  hot-electron  Debve 

oh 

length.  The  oscillation  of  the  potential  barrier,  however,  is  a 
relaxation  oscillation,  the  return  toward  the  target  being  much  more  rapid 


than  the  outward  expansion  (Birdsall  and  Bridges,  1961)  and  thus  has  i 


J 


broadband  frequency  spectrum. 

Using  the  capacltlvely  limited  current  estimate  of  Eq.  (3.61),  with 
potential  barrier  height  |e$  |  ~  3  T  (eV),  this  gives 


„  «.  M  1/2  n.  T.  1/A 

f  (CH.)  -  0.8(i^_ni)  [(  ■  .3)f-ioo^V^l 

r  1 0  cm 


as  a  typical  fundamental  frequency  of  the  broadband  spectrum  of  radiation 
from  the  vibrating  virtual  cathode.  From  Eq.  (3.62)  and  Eq.  (3.9),  this 
should  scale  roughly  as  laser  Intensity,  and  with  wavelength  as 


[Xl/3f(X)]_1  . 


The  power  radiated  by  a  charge  q  oscillating  over  a  distance  d  Is 


-IS  2  2  2  2  A  -A 

PT(W)  -  1.1  x  IQ  q  (Coul  )d  (m  )o>  (s  )  .  (3.1 1A) 


If  the  charge  q  is  the  hot  electron  content  of  the  space-charge  cloud 


"  '  e,co  I  Vn  (x  ■  0) 


(3.115) 


from  Eq.  (3.36),  one  has 


,-7  ,  r0  s2. 


n.  (x  =  0)  T.  1/2 
h  sr  h  ->i 


-  1.88.10-  f_)  l|*J  -(3.116) 

1 0  cm 


For  d  “  1/2  r  ,  l.e. ,  a  50Z  oscillation  in  the  thickness,  r„,  of  the 
q  M 

space-charge  layer,  and  for  U  |  -  e|$  | /T  *  3,  one  has 

m  ra  n 


A .  .  . .  ,  .-  v 


A 

(scaling  as  dz).  Using  the  capacitively  limited  current  to  determine  the 
oscillation  frequency  then  gives  the  power  scaling 


PT(W)  ~ 


6.7  x  io"4  (L“T[ 


nh(x-°)1r 

— ^rJl 

10  cm 


*h 

I  keV 


3  rn  4 

)  (-j— — ) 


(3.118) 


for  50Z  oscillation  in  thickness  of  the  layer,  where  Is  the  rlsetlme  of 

the  laser,  n  (x*0)  is  the  hot  electron  density  of  the  critical  surface, 
n 

and  rQ  is  the  radial  extent  of  the  space-charge  layer,  approximately  equal 
to  the  laser  spot  size.  In  fact,  the  power  radiated  by  the  space-charge 
layer  oscillation  is  generally  somewhat  larger  than  this  because,  scaling 
as  u  ,  It  Is  weighted  toward  the  higher  harmonics  of  the  sawtooth 
oscillation  waveform,  whereas  Eq.  (3.118)  is  calculated  based  on  the 
fundamental  frequency. 

The  radiation  electric  field  corresponding  to  virtual  cathode 
oscillation  is 


E  (V/cm) 
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1/2 

27 .A  [  P_(W) ]  [ 


sin20  i 

“2 - T~\ 

r  (cm  ) 


1/2 


(1.119) 


measured  near  the  ground  plane  at  distances  larger  than  the  wavelength 

c/f  ,  this  is 
osc* 

Eyc(V/cm)  -  9 .47  x  10_2[PT(W) ] 1/2(r/l  ra)_1 


i  *  in“3  (2*5  n«H  nh(x“0)  i1/2/  Th  s3/2,  r0  ?  ,r  -,"1 

)[7nTr;-^  Itsp  W 
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mm'  'TV 

(3.120) 


or,  using  Eq.  (3.9)  for  Tu  as  a  function  of  IX  , 
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(3.121) 


3/4 

which  scales  as  1^  /f(X)  if  n^  Is  scaled  according  to  Eq.  (3.62)  and  If 

Th  “  dLx2)1/3‘ 

This  is  to  be  compared  with  the  quaslstatlc  field  due  to  the 
unneutralized  "beam"  space  charge  of  the  longer  plume  of  escaping  fast 
electrons  (Eq.  3.107) 


,  I  (A)[L(cm)]‘ 

E  .  (V/ cm)  ~  2.4  x  10  -2 - J  . 

plume  f 


{e<J>  (keV)] 
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r-r3 

m' 


(3.122) 


With  ~  1A,  e^m  ~  1  keV,  and  a  20  cm  plume,  this  lower -frequency  field 

at  the  detector  (r  *  1  m)  is  of  order  0.1  V/ cm,  with  rise  time  of  order 

3  ns.  Assuming  e$  «  T,  *  (I  X^)  * ^  from  Eq.  (3.9),  *  x2e4>  (L  -  vt), 

ran  m 

3/2 

and  current  1^  scaling  as  1^  ,  the  quaslstatlc  plume  electric  field  would 

scale  as 


E  i  * 

plume 


t2(Il)5/SX1« 


(3.123) 


with  t  the  laser  pulselength,  whereas  the  radiation  field  from  virtual- 

3  1/2 

cathode  oscillations  would  appear  to  scale  as  (nhTh^  '  /Tr ’  *,e*  »  as 
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E 

VC 

with  the  risetime  of  the  laser  pulse  and  with  f(X)  presumably  only 
weakly  dependent  on  X .  Based  on  the  observed  plume  current  it 

would  appear  that  the  fundamental  GHz  field  from  vc  oscillation  has 
smaller  magnitude  (~  .5  mV/cra  at  laser  intensities  I  ~  10  W/cm  )  than 
the  monopole  field  of  the  plume  unless  the  effective  plume  length  is 
~2  cm  or  smaller. 

Regardless  of  the  scaling  assumed  for  n^,  Th»  and  1^,  one  should 

have 


E 


plume 


(3.12A) 


and 


vc 


(3.125) 


As  discussed  in  Chapte4  A,  Faraday  cup  arrival -time  data  show  Th 

1/3 

scaling  even  more  weakly  with  IL  than  1^  ,  and  magnetic  loop  signals 

~1  75 

appear  (if  functioning  correctly)  to  show  1^  *  »  while  the  electric 

field  data,  somewhat  more  ambiguous,  could  be  interpreted  as  either 
E  *  IL  or  E  «  IL  . 


Chapter  4 


DATA  ANALYSIS 

Of  the  1200  shots  taken  during  the  experiment,  a  highly  selective 
subset  was  chosen  for  detailed  analysis  through  mutual  discussion  between 
AFWAL  and  PRC.  The  shots  represented  parametric  scans  in  intensity  at 
le-6  torr  and  at  3e-3  torr  with  a  full  data  compliment,  and  pressure  from 
le-6  torr  to  l  torr.  The  shots  chosen  had  a  full  diagnostic  set  and  a 
nominal  laser  pulse  waveform  as  shown  in  Figure  4.1.  Data  from  these 
shots  were  subsequently  reduced  at  AFWAL  by  Technology  Scientific  Services 
and  the  data  were  differentiated  and  integrated  per  PRC  request.  Whenever 
possible,  five  data  sets  were  chosen  for  each  parametric  point. 

The  graphical  data  were  then  tabulated  for  subsequent  plotting  as  a 
function  of  intensity,  pressure  or  other  dependent  variable  to  yield  the 
scaling  information  for  comparison  with  the  theoretical  predictions. 

Where  sensible,  least  squares  linear  fitting  was  performed  on  the  data 
sets  to  derive  ei^>lrical  scaling  relationships.  Other  manipulations  of 
the  data  are  described  on  a  case-by-case  basis. 


4.1  GENERAL  WAVEFORM  ANALYSIS 

The  observed  waveforms  for  the  Intensity  data  are  shown  In  Figure 
4.2.  Data  from  the  Electric  Field  Detector  (d/dt  of  the  Field)  DD2  are 
shown  at  two  pressures,  le-6  torr  and  3e-3  torr,  as  a  function  of  the 
attenuator  thickness  (and  thus  Intensity).  The  transmissions  of  the 
attenuators  are  given  by: 
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Thickness 


Transmission 


0 

mm 

100% 

1 

mm 

64% 

3 

mm 

32% 

6 

mm 

9% 

The  transmission  of  the  salt  vacuum  window  was  92X. 

The  individual  photos  also  show  the  volts  and  time  per  division  of 
the  oscilloscope  and  the  shot  number  of  the  data.  Except  for  shot  651, 
which  is  the  only  low  intensity  le-6  torr  shot  with  a  full  data  set,  the 
other  shots  are  representative  of  the  data  used  in  the  analysis. 

The  DD2  data  at  le-6  torr  and  low  laser  intensity  are  very  noisy 
and  a  reproducible  pattern  is  hardly  discernible.  As  intensity  is 
increased,  a  negative-going  signal  is  observed  that  increases  with 
increasing  intensity,  yet  retains  high  frequency  components.  This  high 
frequency  ringing  is  from  plasma  frequency  oscillations  in  the  expanding 
blowoff  plasma  and  probably  represents  the  low  frequency  end  of  a  spectrum 
of  oscillations  that  generally  exceeded  the  bandwidth  of  the  recording 
system.  The  noise  in  the  system  made  hand  digit laatlon  less  accurate  than 
would  be  deslreable  for  deriving  scaling  information  and,  in  retrospect, 
it  would  have  been  better  to  passively  Integrate  the  d/dt  signals  to 
obtain  the  fundamental  field  parameters. 

At  3e-3  torr,  the  DD2  signal  Is  stronger,  and  the  high  frequency 
oscillations  have  been  damped  by  Intervening  plasma  In  the  vacuum  tank. 
Again  the  peak  signal  values  Increase  with  Intensity,  and  the  duration  of 
the  main  signal  Is  the  same  as  the  laser  pulse  length.  There  are 
differences  In  the  fluctuations  for  shots  974  and  9S8  that  are  not 
analysed  since  they  occur  late  In  time.  hate  oscillations  are  observable, 


but  are  not  considered  in  the  data  analysis. 

At  both  pressures  the  data  traces  froa  Faraday  Cup  #1  (FC1) 
increase  with  Intensity  and  are  slailar  in  shape.  The  initial  bump  In  the 
traces  is  due  to  the  differentiated  pyroelectric  signal  that  was  used  as  a 
timing  aarker.  The  electrons  arrive  at  the  detector  about  20  ns  after  the 
leading  edge  of  the  laser  pulse.  There  was  a  3.4  ns  cable  delay  that  oust 
be  subtracted  from  the  timing  shown  for  FC1  to  compute  the  electron 
energies.  The  energies  of  the  first  observed  electrons  are  of  order  1  keV 
at  the  FC1  station  30  cm  from  the  target.  The  FC1  waveforms  at  3e-3  torr 

have  a  slower  decay  than  the  le-6  traces. 

The  PMM  traces  measured  the  time  derivative  of  the  current  through 
a  loop  at  a  radius  of  3  cm  from  the  target  In  the  ground  plane.  The 
traces  at  hard  vacuum  are  of  the  same  duration  as  the  laser  pulse,  while 
those  at  3e-3  torr  show  a  significant  tall  that  Is  longer  than  the  laser 
pulse  by  a  factor  of  3  or  more.  PRC  believes  this  Is  due  to  flashing  of 
the  detector  or  In  the  ground  plane.  The  detector  had  an  Inductance  of  2 
nH  and  the  L/R  decay  time  of  the  trace  leads  to  a  ’-esistance  of  .2  ohms 
which  Is  consistent  with  plasma  resistivities  in  parallel  with  the  SO  ohm 
cable  Impedance.  The  flashing  phenomenon  is  clearly  Indicated  in  Figure 
4.3. 

At  a  pressure  between  le-5  and  le-4  torr  there  Is  an  abrupt  change 

In  the  FMM  waveform  to  the  long  decay  shape.  The  change  occurs  at  the 

peak  In  the  FMM  signal  and  usually  Is  associated  with  a  glitch  In  the 
waveform.  Once  the  waveform  changes  shape  they  are  similar  up  to  the 
highest  pressures  of  .1  to  1  torr,  where  the  waveform  Is  further  modified 
by  background  pressure.  The  f)R2  and  B2  waveforms  Indicate  an  electro¬ 
magnetic  wave  is  observed,  since  the  waveforms  are  similar  at  all 


pressures  and  there  is  no  observable  phase  change  between  the  two.  The 
high  frequency  noise  In  DD2  and  B2  at  le-6  torr  Is  rapidly  damped  at 


higher  preasures  and  the  strength  of  the  wave  Increases  nodestly  at  higher 
pressures.  Long-term  ringing  of  the  ground  plane  is  apparent  after  the 
Initial  lapulse  produced  by  the  laser,  so  data  after  10-15  ns  Is 
contaminated  with  reflections. 

4.2  DATA  TABULATION 

i 

t 

t 

To  obtain  more  quantitative  data,  the  waveforms  and  moments  were  j 

i 

tabulated  at  given  times,  the  tabulated  data  were  then  plotted  against  j 

t 

intensity  or  other  Independent  variable  to  derive  scaling  Inf ormat ion.  \ 

i 

The  points  where  the  Information  was  taken  are  shown  In  Figures  4.4  J 

I 

a,b,c,d.  Figure  4.4a  Is  laser  Intensity  data,  4.4b  Is  DD2  data,  4.4c  Is  | 

FMM  data,  and  4.4d  Is  Faraday  cup  data.  The  Intensity  tabulation  Is  shown 
In  Table  4.1  and  the  pressure  data  are  shown  In  Table  4.2. 

t 

Figure  4.4  shown  the  Intensity  (pyrometer)  data  for  shot  612,  both  | 

as  original  digitized  data  and  also  Integrated  and  differentiated,  as  a 
function  of  time  In  nanoseconds.  This  data  was  tabulated  for  all  shots 
for  comparison  to  the  calorimeter  data  (cal)  and  used  to  derive  peak 
Intensity.  The  peak  pyro  d/dt  was  also  tabulated  for  comparison  to  rate- 
dependent  data,  but  the  data  have  significant  spikes  from  the  digitization 
process  that  need  to  be  smoothed  before  believable  Inferences  could  be 
made. 

Figure  4.4b  shows  the  DD2  data  (t»,  or  time  derivative  of  electric- 
field).  The  original  data  are  spiky  by  nature  and  the  peaks  are  recorded 
for  plotting.  Hore  believable  data  occur  after  integration  to  obtain  the 
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electric  fields,  and  the  peak  electric  field  was  tabulated  at  1DD2.  The 


"action  Integral"  of  the  DD2  signal,  l.e.,  the  time-integral  of 

2 

(dE/dt)  ,  was  tabulated  at  the  peak  field  and  was  used  as  a  consistency 
teat  of  the  digitisation  since  this  should  scale  as  the  square  of  the 
fields  and  is  proportional  to  the  energy  in  the  electric  field. 

Figure  4.4c  shows  the  tabulation  points  for  the  FMM  (magnetic  field 
time  derivative)  data.  The  current  and  the  charge  ( IFMM  and  1IFMM)  were 
taken  at  40  ns  to  avoid  late-time  drifts  in  the  detector  that  occurred 
after  the  laser  pulse  was  off. 

Figure  4.4d  shows  the  tabulation  points  for  the  FC1  data.  The  zero 
time  is  taken  as  the  peak  of  the  rlsetime  of  the  fiducial,  which  is  a 
differentiated  pyro  signal.  This  places  t  -  0  at  the  leading  edge  of  the 
laser  pulse.  Subsequent  times  are  taken  from  the  differentiated  waveform. 
Tfcl  la  where  the  trace  goes  through  zero  and  represents  the  leading  edge 
of  the  observable  signal  (fastest  observed  electrons).  TfcIIn  is  the 
Inflection  time,  i.e.  ,  when  the  rate  of  rise  of  the  FCl  signal  Is  maximum. 
TFC1M  is  the  time  of  the  maximum  Faraday  cup  signal.  All  the  signal 
timings  Include  the  measured  TDR  times  for  cable  lengths. 

Table  4.1  shows  the  data  taken  at  both  l  .e-5  and  3.e-3  torr  as  a 
function  of  the  filter  thickness.  The  laser  data  was  then  used  to  derive 
the  peak  intensity  using  the  measured  spot  st  -.e  and  spatial  distribution 
of  the  energy.  The  intensity  of  the  peak  was  computed  using  a  gausslan 
approximation  to  the  data  which  leads  to  a  factor  of  1.93  between  the  peak 
and  the  average  Intensity.  Thus  If  the  reader  needs  the  spatially 
averaged  Intensity,  he  can  divide  by  1.93  (or  subtract  .29  from  the  log) 
of  the  tabulated  data.  Where  data  Is  missing,  it  was  either  not  suitable 
for  reliable  reduction  or  was  missing. 
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Table  4.1 


Intensity  Tabulation 


allot  pr 
#  torr 


Pt  PT  t  Xpy  py/at  lot  an 
m  ns  T*na  Mv/a*W/an2 

•16  354  10.6  3.8  300  1?# 

.29  410  10.4  4.4  520 

‘  2  11. - 

0  10. 

3  10. 


i/a 


IfH  life  dd2  Id2  Ud  fan 
in  pa  uooul  ref  T  4.2  nf 


t fan  tfo 

ns  ns 


tfol  Ifo 
ns  naoul 


64  10.1 
44  10.2 
0  10.5 
0  10.6 
0  10.4 


1% 

500 


is 


0.099 

0.111 

0.086 


1 

1  0.079 
4  0.086 


2.7 

2.1 


0.04  — 
0.05  260 


17 


132  366 

182  530 
1.7  0.05  170  147  362  200 
1 .7  0.06  230  1  80  670  210 
1.7  0.04  1  80  1  50  420 


30.0  22.0  28.5  1.52 
26.0  20.0  24.0  2.50 


29.0  20.5  24. 
28.0  20.0  “ 


24.5  2.01 
24.0  2.50 


29.0  21.0  24.5  2.15 


0.9  0.02 
1  .0  0.02 
0.8  0.02 
0.8  0.02 
0.8  0.02 


17 


40  1 


<j  I?  9*  1 

n  ini ! 

56 


27 


1.1-06  1.38 

400 

10-3 

4.2 

440 

1.22 

0.015 

0.2 

0.01 

I  68 

71 

3"  1  •« 

m 

11.0 

4.7 

700 

1.18 

0.015 

0.3 

0.01 

1.29 

10.4 

4.5 

595 

1  .18 

0.021 

0.2 

0.00 

4  28 

17 

i$£°6  1,37 

420 

10.6 

4.4 

556 

0.69 

0.017 

0.0 

0.78 

-  - 

- 

0.0  20.0 
2.0  20. 

2.0  21 . 

2.0  22.0  25. 
3.0  22.0  26. 


.0 

.0 

.0 

5 

0 


1.47 
1.71 
1.68 
1.43 
1  -57 


2.0  20.0  27.0  1.0| 


0  0.99 
33-5  23.0  25.5  0.43 


IO.3  4.5 
10.8  3.6 

11.2  4.3 
3  10.4  3.f 


10.8  4.0 

10.2  3.6 

11.2 
11.0 
10.2 


405  3.8-03 

401  6 an 

409 

410 
413 


1.06  358  9.6 

!:!l  iio ,0* 

!:« 131  !o. 


3.700  218.0  4.80  120  244  775  172 

1.63  2.900  278.0  5.20  27  218  165  156 

1.66  3.OOO  228.0  4.80  60  191  121 

1.68  2.000  135.0  2.70  60  190  61 8  84 


1.46  2.500  159.0 

1.47  2.500  156.0 


3:8  ?  Il^i. 


1.43  3.OOO  162.0  4.O0  13  117  282  84 

1.45  2.400  144.0  2.70  6  102  232  90 

1.47  2.300  130.0  2.40  8  110  253  99 

0.90  - 

I.63  _ 

?:»  = 


11.0  4.; 


5.0  22.5 

~.o  22.5  „ 

.0  22.0  26.0  2.31 

2.0  20.0  29.5  1 .57 


.0  2lfo 


.0  22.0  29.5  1 .48 

0.0  1.14 


.0  22.5 
.0  22.3 


.  .5  22  .. 
4.0  21  .5 


2.0  1.4( 
2.0  1.' 
1.0  l.i 


1.20  0 


1.13  0. 
1.18  0. 


1.21  1. 
1.09  0. 


0.0  0.96 
0.0  0.95 
.0  1.22 
.0  1  .05 


2.0  22.0  24.0  1.17 


0.131 


0.6’ 

0.65  0.117 
0.66  0.114 

0.60  O.I35 
0.59  0.125 


7.4  O.13 
1.10 


?:?  °- 


.0  0.10 
7.2  0.12 
7.7  0.13 


-  41  29  27 

1  25  29  27 

:  S  H  a 

-  53  54  24 


0.5  22.0  24.0  0.42 
2.0  23.5  25.0  0.42 
2.0  23.0  25.0  0.45 
4.0  23.0  28.0  O.36 
3.0  24.0  26.0  0.31 
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Table  4.2 


Preaaure  Tabulation 


preaa  cal  py  pk  py  t  Ipy  r 
torr  J  wr  ns  v*na  f 

1.8-06  1.25  410  10.3  4.40 

1.23  417  10.8  4.50 


Ifaa  liras  dd2  Idd2  Udd2 
Aapa  pcoul  7/a/a  V/a  p?* 2s/a“2 


1 .8-05  1 .20 

1 .42 
1.38 


1 .8- 04  1 .52 

1 .09 
0.91 
1 .03 
0.97 

3.8- 04  1.18 

1 .10 


1 .8-03  1 .05 

0.90 
0.93 
0.98 


3.8-03 


6.8-03 


1.8-02  1.20 
1.23 
1.18 


2.8-02  1.04 

1.13 

1.14 


6.8-02  1.02 
1  .09 
1  .01 
1.16 
1  .00 

1.8-01  1.21 


1.8*00  1.15 

1.15 


:!  P 

.3  4.0C 


^  1  0  •  0  J  «  -TV 

}0  10.0  3.50 
>4  10.4  3.80 


417  10.8 
370  10.1 
382  1  0.1 
387  10.3 


83  10.1  4.10  310 
16  9.8  5.10  360 

40  7.8  4.70  320 

158  10.5  4.00  185 

70  10.5  3.90  200 


487  9.2  4.80 
360  10.8  3.90 
300  10.0  2.90 

339  10.9  3.40 
328  10.3  3.25 

414  9.8  3.90 

!8:S  \:H 

340  10.1  3.3O 
365  10.0  3.70 

330  10.7  3.50 
310  10.2  3.IO 
300  10.5  3.10 
328  10.4  3.40 

310  11.0  3.30 

405  10.8  4.10 
450  10.3  4.40 
410  10.8  4.40 
365  10.4  3.7O 
360  10.4  3.70 

420  10-3  4.20 
460  10.4  4.70 
438  10.2  4.40 
505  9-7  4.80 
422  10.8  4.50 

405  11.2  4.30 

a  10. 4  4.20 
10.0  3.90 

10*.!  3*.W 

326  10.4  3.40 
340  10.9  3.70 
3?4  10.6  3.40 


10  0.19 
60  O.39 
20  0.21 


1.44  53.OO  I.36 

0.75  45.00  1.00 

0.69  39.00  1.00 

1.04  72.00  I.95 
0.72  43.00  1  .20 

2.40  135.00  3.9O 
1.98  126.00  3.45 

2.00  119.00  3.40 

2.90  163-00  4.70 

3.00  171.00  4.80 

5.90  388.00  11.00 

5.00  350.00  10.00 
4.50  300.00  8.00 

4.50  394.OO  8.20 

4.00  380.00  8.00 

1.40  114.00  6.00 

.90  240.00  7.00 

.70  23Q.OO  7 -00 
.30  208.00  5.50 

.10  107.00  6.00 


0.07 

105 

172 

360.0 

0.05 

— 

— 

— 

0.04 

144 

131 

90 .0 

260.0 

0.04 

137 

165 

0.13 

— 

— 

0.08 

164 

272 

300.0 

0.12 

0.07 

225 

240 

1^9 

760.0 

760.0 

0.09 

145 

300 

450.0 

0.07 

181 

338 

630.0 

1.36 

400 

512 

BOO  .0 

1  .00 

225 

310 

360.0 

1  .00 

143 

182 

120.0 

1.25 

135 

158 

80.0 

5.00  290.00  9.00 

6.6  0  350.00  10.00 
6.00  320.00  9.00 

7.00  338.OO  10.00 
5.40  324.00  9-00 


9.0° 

16.OO 


6.15  447.00  12.00 
6.80  452.00  n.OO 
6.70  443.00  13.00 


192  213 

m  & 

124  114 

104  102 


'll  108 

179  175 


4.60  371.00  10.00 
4.90  410.00  11 .00 

5.60  49O.OO  13.00 
5.45  420.00  12.00 
6.10  477.00  13.00 


210  l89 

180  156 

185  221 

180  158 

200  114 

111  il! 


200  118 
if.  15 
111  ,« 


\u  ,n 

110  175 

IS  8 


317  10.2  3.20  168  4.00  479.00  10.00 

325  10-I  3-5°  359  4.35  480.00  11.00 

,05  11.4  j.M  10$  4.50  543.00  12.00 

20  11.2  .60  164  4.40  527.00  12.00 

336  10.0  3.3O  156  4.15  489.00  12.00 

378  10.7  *.20  190  3-80  516.00  11.00 

810.8  3.70  220  3-50  490.00  10.00 

18:11:8  118  118:88  »:« 

10.8  3.80  200  3.06  495.00  34.00 

357  10.8  3.90  1  30  0.37  90.OO  4.00 

350  10.3  3.00  1*3  0.38  87.00  10.00 


190  3-80  516.00  11 .00 

220  3-50  490.00  10.00 


270.0 

400.0 

90.0 

140.0 
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4.3  ANALYSIS  OF  TABULATION  -  INTENSITY  SCALING 

The  calorimeter  data  was  plotted  against  the  integral  of  the 
pyroelectric  detector  to  check  linearity  of  the  data,  since  the  intensity 
was  computed  using  the  calorimeter  energy  and  the  pyro  FWHM.  The  results 
are  shown  in  Figure  4.5  and  demonstrate  that  the  system  was  well  behaved. 

A  linear  least-squares  fit  goes  nearly  through  the  origin  and  leads  to 
confidence  in  the  intensity  numbers.  Figure  4.6  shows  the  calorimeter  vs 
the  pyro  d/dt  signal.  The  curve  is  nearly  flat  since  the  energy  is 
proportional  to  the  amplitude  of  the  laser  rather  than  to  an  increase  in 
the  pulse  length.  Thus  the  basis  for  "subsequent  characterization  of 
laser  Intensity  by  calorimeter  data  in  this  study"  is  well  founded 
statistical ly. 

The  DD2  data  are  plotted  in  Figure  4.7  in  both  semilog  and  log-log 
forms.  The  semilog  plot  shows  a  threshold  behavior  at  30  Gw/cm^,  but  when 
plotted  log-log,  the  evidence  for  power  law  behavior  is  more  apparent, 
the  power-law  dependence  is  consistent  with  the  current  data  that  will  be 
shown  subsequently.  The  data  for  DD2  are  noisy  and  if  the  lowest  point  in 
the  le-6  data  is  ignored,  there  is  a  substantial  discrepancy  between  the 
behavior  of  the  le-6  and  3e-3  torr  data.  The  data  at  3e-3  torr  scale 
roughly  as  the  1/2  power  of  intensity;  however,  if  the  lowest -intensity 

data  are  ignored,  the  le-6  torr  data  scale  as  the  3/2  power.  The  1/2 
power  scaling  is  roughly  consistent  with  the  theoretical  modeling,  and  the 
experimental  data  at  3e-3  torr  is  strongly  for  the  weaker  scaling.  The 
lowest-intensity  point  at  le-6  torr,  if  valid,  would  discount  the  3/2- 
power  scaling  at  le-6  torr  as  well.  This  is  a  major  discrepancy  that  we 
are  unable  to  resolve. 
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In  Figure  4.8  the  Integrated  DD2 ,  i.e. ,  the  electric  field  data, 
are  shown.  This  integration  presumably  smooths  some  of  the  numerical  and 
detector  noise,  but  the  difference  still  exists  between  the  le-6  and  3e-3 
data  scalings.  The  combined  data  were  fitted  with  a  least-squares 
approximation  and  the  data  scale  as  the  1/2  power  of  intensity. 

The  FMM  data  and  its  moments  are  plotted  in  Figure  4.9  at  le-6 
torr.  The  data  are  well  behaved  and  clearly  power-law  in  nature.  The 
current  and  charge  are  scaling  as  intensity  to  the  1.5-2  power.  Since  E 
and  B  (i.e.  ,  current)  should  scale  together,  this  lends  credence  to  the 
vacuum  case  scaling  about  to  the  3/2  power.  But  as  discussed  in  section 
3.4  (and  subsequently  in  section  5.1)  such  a  scaling  does  not  appear 
consistent  theoretically  with  either  the  weak  dependence  of  hot  electron 
energies  or  with  1/2-power  scaling  of  the  electric  fields.  The  FMM  data 
at  3e-3  torr,  shown  in  Figure  4.10,  also  scale  as  intensity  to  the  3/2 
power,  in  spite  of  the  concerns  about  the  detector  flashing  behavior 
mentioned  before.  This  is  indicative  of  the  detector  flashing  at  the  time 
of  the  large  trailing  edge  dl/dt,  thus  preserving  the  largest  amplitude  of 
FMM  before  L/R  decay. 

The  Faraday  Cup  data  are  plotted  as  a  function  of  peak  intensity  in 
Figure  4.1.  Since  there  was  little  difference  in  FC1  as  a  function  of 
pressure,  data  for  both  pressures  are  plotted  together.  The  energies  of 
arriving  electrons  scale  as  the  .1  power  of  intensity,  while  theory  would 
predict  a  .2-. 4  dependence.  Given  the  scatter  in  the  data,  this  is  not 
far  from  the  observations.  At  30  Gw/cm^,  the  mean  energy  for  the  leading 
edge  is  1.6  keV,  the  mean  inflection  energy  is  1  keV  and  the  mean  energy 
at  the  peak  is  630  eV. 

The  FC1  currents- and  charge  were  plotted,  Figure  4.12,  for  both 


pressures,  showing  that  the  FC1  current  scales  to  the  3/4  power  of  the 
Intensity,  which  In  turn  Implies  a  1/2  power  scaling  with  the  FMM  since 
FMM  scales  as  the  3/2  power  of  Intensity.  This  relation  is  confirmed  in 
Figure  4.12.  This  Is  convincing  experimental  evidence  for  beam-like 
propagation  that  Is  intensity-dependent,  rather  than  a  uniform  cosine 
distribution  of  the  charge.  There  were  no  Faraday  Cups  at  other  azimuthal 
angles,  so  we  do  not  have  direct  experimental  evidence  of  the  angular 
distribution,  but  must  rely  on  the  inference  from  this  plot. 

4.4  ANALYSIS  OF  TABULATION  -  PRESSURE  SCALING 


Since  the  Integral  of  DD2  is  the  Electric  Field  and  a  fundamental 
quantity  as  well  as  being  smoother  than  the  DD2  data,  the  pressure 
relations  were  plotted  In  field  and  FMM  current  ( 1FMM) .  Figure  4.14  is 
the  electric  field  and  current  scaling  as  a  function  of  pressure.  In  this 
plot  all  data  points  in  Table  4.2  are  shown.  There  is  considerable 
scatter  In  the  semi  log  plot,  with  f actor-of -two  fluctuations  common. 

Since  there  were  fluctuations  in  the  laser  intensity  that  were  small,  the 
data  were  linearized  to  the  laser  calorimeter  energy  and  averaged  before 
replotting  as  In  Figure  4.15.  This  figure  also  shows  the  regions  where 
Ion  and  thermocouple  gauges  were  used  to  measure  the  pressure.  On  sample 
wMsurfwnts  In  the  overlap  region  these  devices  Indicated  a  discrepancy 
>f  70-fold;  this  region  Is  approximately  indicated.  Thus,  the  data  points 
In  and  around  this  region  may  be  translated  significantly  from  their 

apparent  position.  For  example,  the  points  at  le-3  torr  could  be  at  Ie-2 
torr  and  then  the  curves  would  be  continuous  rather  than  showing  the  bump 
behavl or . 
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The  electric  field  data  show  a  graduate  rise,  betwwen  le-6  and  le-5 
torr,  of  a  factor  of  two.  Although  this  appears  statistically 
significant,  it  may  he  due  to  improved  digitization  of  the  original 
waveforms,  since  the  very  high  frequency  ringing  is  damped  rapidly  ns  the 
pressure  rises.  There  is  no  ready  theoretical  explanation  for  this  rise 
in  E-field  at  such  low  pressures  and  an  argument  can  be  thus  made  that  the 
electric  field  scaling  should  be  flat  until  le-3  to  le-2  torr  where  the 
effects  of  electron  dE/dx  begin  to  take  hold  in  the  background  gas.  After 
that  the  electric  field  intensity  drops  nearly  exponentially  with 
pressure,  which  is  the  expected  classical  behavior  for  electrons  in  a 
neutral  background  gas.  The  dE/dx  of  the  electrons  scales  as  1/E;  the 
ranges  vary  as  a  function  of  the  energy  in  the  distribution  at  least  until 
the  peak  of  the  Ionization  cross  section  at  100  eV,  Figure  4.16.  The 
lowest  energy  of  electrons  observed  in  the  experiments  is  300  eV,  although 
electrons  must  exist  to  eV  energies.  The  range  of  the  fast  electrons,  as 
discussed  in  section  3.7,  is: 


e  folds  at 


Energy 

Range 

le-2  torr 

1  keV 

7  cm/P(torr) 

.09 

.3  keV 

2  cm/P(torr) 

.30 

.1  keV 

.7  cra/P(torr) 

.9 

Since  electrons  are  present  down  to  relatively  low  energies,  the 
low  energy  electrons  from  the  distribution  are  slowed  first  as  the 
pressure  rises,  and  this  has  the  effect  of  reducing  the  field.  As  the 
pressure  rises  to  .1  torr,  the  range  of  even  the  faster  electrons  becomes 
comparable  to  the  system  size  and  these  electrons  are  also  slowed,  further 
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reducing  the  field  Intensity 


If  the  Paradsv  Cups  had  Seen  present  for 


Che  pressure  scan,  these  results  could  have  Seen  verified  with  direct 
electron  energy  atasurewnti. 

Figure  4.17  shove  the  pressure  dsta  on  a  lop-log  plot  to  highlight 
the  current  dependence  on  the  pressure.  The  current  Increases  until  .1 
torr  when  the  range  of  the  electrons  becomes  sea  1 1  coepared  to  the  evetee 
size  whereupon  the  signal  decays  towards  zero,  lore  quantitative 
statements  are  not  possible  since  the  FMM  detector  showed  evidence  of 
flashing  above  3e-4  torr. 


4.5  EXPERIMENTAL  SCALING  SUMMARY 


Analysis  of  the  experimental  data  leads  to  scaling  rules  for  the 
currents  as  le-6  torr  and  for  the  Faraday  cup  currents  at  both  pressures. 

FMM  current  ~  intenslty^^ 

FC1  current  ~  lntenslty^^ 

FC1  energy  ~  Intensity’* 

FC1  current  ~  FMM  current*^ 


For  the  electric  fields  there  is  conflicting  data.  At  le-6  torr 
the  bulk  of  the  data  appears  to  scale  as 


Field  ~  Intensity 


3/2 


3/2 


which  is  consistent  with  theory  if  the  current  also  scales  with  1^  ; 
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hoiMvtr ,  ow  data  point  at  the  vary  lowest  intensity  could  be  construed  to 
bring  the  scaling  to  the  1/2  power,  consistent  with  the  3e-3  torr  data 
which  caper iaent al ly  acalea  aa 

,  .  .  .  1/2 

field  ~  Intensity 

This  Is  roughly  consistent  with  theory  If  the  current  does  not  scale  aa 
1^^  but  as  some  considerably  weaker  power  of  1^,  e.g. ,  1^^  as  discussed 
at  the  end  of  Sec.  3.4.  The  electric  field  observations  were  eade 
difficult  by  high  frequency  components  at  low  pressure,  yet  the  cleaner 
data  does  not  agree  with  the  theoretical  expectation  that  the  field  should 
scale  nearly  the  same  as  the  current. 

The  scaling  of  field  with  pressure  Is  well  explained  by  classical 
dE/dx  of  the  electrons  in  the  background  gas  leading  to  reduced  charge 
separations.  The  electric  field  Is  reduced  20Z  per  decade  of  pressure 
between  le-5  and  1  torr.  This  Is  an  exponential  reduction  In  the  field 
strength  with  pressure  and  consistent  with  the  dE/dx  calculation  In  the 
text.  This  effect  coupled  with  beam  propagation  calculations  are  well 
founded  theoretically,  and  experimental  Inferences  are  consistent  with 
this  model.  Direct  observation  of  the  Isotropy  of  the  electrons  and  the 
change  of  energy  with  pressure  are  not  available,  since  the  Faraday  cups 
were  not  fielded  on  the  pressure  scans  or  at  varying  azimuthal  angles.  In 
future  work,  placement  of  identical  cups  at  varying  angles  and  maintaining 
these  cups  throughout  the  experiments  is  recommended.  In  addition,  DD2 
and  FMM  data  should  be  passively  electronically  integrated  before 
recording,  In  order  to  eliminate  digitization  difficulties. 
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Figure  4.1 
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THE  ELECTRIC  FIELD  SCALING  FOR  POINTS  AT  IE-6  AND  3E-3  TORB 
PLOTTED  TOGETHER  AS  A  FUNCTION  OF  PEAK  INTENSIIV.  THE 


LOG  PEAK  INTENSIIV  (  GH/cnfi) 


LOG  PEAK  INTENSITY  (  GH/cm2) 

0f  mm*100  Alf  mh#100  01  IftiM*le4 


ALL  FMM  MOMENTS  SCALE  AS  INTENSITY  TO  THE  3/2  POWER 
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Chapter  5 


SCALING  OF  OBSERVABLES  WITH  LASER  INTENSITY,  WAVELENGTH, 
AND  BACKGROUND  GAS  PRESSURE 


5.1  INTENSITY  SCALING 


5.1.1  Total  Current 

The  fast-electron  current  1^,  as  measured  by  a  magnetic  loop  around 
the  target  In  the  target  plane,  scales  with  the  laser  intensity,  1^,  as 


I 


b 


.3/2-2 


(5.1) 


-3 

with  the  3/2  observed  at  P  •  3  *  10  torr  and  the  2  observed  at 
P  ■  10  6.  (See  Figs.  4.9  and  4.10.)  There  appear  to  be  no  theoretical 
calculations  in  the  literature  giving  any  scaling  laws  for  the  rate  of 
fast-electron  production  as  a  function  of  laser  intensity,  although  there 
are  many  giving  scaling  laws  for  the  fast-electron  temperature. 


nnnminmnunf«iniiniKnvtu«uininviMW 


5.1.2  Current  Density  at  the  Faraday  Cup 

The  current  density,  measured  by  Faraday  cup  10°  from  the  laser 
axis  at  30  cm  from  the  target,  scales  approximately  as 

J(z  -  29.54,  r  -  5.21)  «  lj/4  (5.2) 

—6  -3 

at  both  pressures  where  the  scaling  was  evaluated  (10  and  3  *  10 
torr),  indicating  that  the  effective  area  over  which  the  beam  or  plume  of 
fast  electrons  was  spread  varied  as  the  square  root  of  the  current: 

J  «  I*/2  *  Agff  ( z  -  29.54)  -  lj/2  (5.3) 

(see  Figures  4.12  and  4.13). 

If  the  emission  always  occurred  with  a  cosine  angular  distribution 
as  is  sometimes  postulated,  one  would  expect  J  *  I, .  Thus,  there  is  some 

D 

current-dependent  beam  spreading  of  the  fast-electron  plume,  and  the  beam 
of  fast  electrons  is  narrower  than  a  cos8  distribution  at  least  at  the 
lower  currents.  The  actual  beam  radius  at  the  detector  depends  on  both 
the  initial  angle  spread  48^  and  on  the  space-charge  spreading.  For 
A8q  small,  one  can  fit  the  results  of  Table  3.1  to  a  form 

a2(cra2 )  -  200  [ Ib(A)/E(keV)3/2 ) 1  * 1 1 7  (5.4) 

at  the  detector  (a  is  the  beam  radius).  Since  this  value  of  a  is 
typically  larger  than  the  5.2  cm  detector  radial  position  (measured  from 
the  axis)  by  only  a  factor  of  2,  and  since  the  plume  probably  has  a  bell¬ 
shaped  distribution  of  current  density  with  width  ~a,  the  actual  current 


density  Measured  at  the  Faraday  cup  could  scale  In  a  complex  way,  e.g., 


J  -  *-1Iwa‘2[l  +  (5.21/a)2)'2 

D 


(5.5) 


for  a  Bennett  distribution  with  radius  a  given  approximately  by 


a2  ~  a2(A0 )  +  200  { I  /E3/2 ) 1  -1 1 7 

0  b 


(5.6) 


(with  a  In  cm,  In  A,  and  E  In  keV).  When  the  two  terms  in  Eq.  (5.6) 


are  comparable,  J  can  scale  roughly  as  I*'  over  an  order  of  magnitude  in 

b 


5.1.3  Electric  Field 


The  peak  electric  fields  seen  by  the  "D  dot"  detector  1  m  from  the 


source  near  the  ground  plane  (z  -  2  cm)  scale  with  laser  intensity  as 


1/2  -3 

1^'  at  P  -  3  *  10  torr  (Fig.  4.3),  while  in  vacuum  the  dependence  on  1^ 


is  less  clear  but  could  be  stronger,  e.g., 


e  «  i; 


(5.7) 


Theoretically  the  low  frequency  part  of  E  should  scale  as  1^  until 


space-charge  neutralization  by  ionization  sets  in  at  ~  10  torr  (actual 


pressure  depending  on  It,)*  The  higher  frequency  part  of  E,  which  is  seen 


from  the  data  traces  (Figures  4.2  and  4.3)  to  dominate  at  low  P,  should 


scale  as  fluctuation  amplitudes,  e.g.. 


V  *  E 


Si  /t 
net  osc 


(5.8) 


with  t  *(osc)  * 

-3 

At  P  -  3  *  10  torr,  these  fast  oscillations  are  somewhat 

suppressed  and  the  E  field  (Integral  of  0D2  in  Figures  4.1  and  4.2)  is 

dominated  by  the  low-frequency  component,  which  is  presumably  proportional 

1/2 

to  the  net  current  Inet.  Nonetheless,  the  observed  1^  scaling  of  the 

1/4 

electric  field  is  consistent  with  theory  only  if  I  *  I.  at  this 

net  b 

pressure. 

5.2  WAVELENGTH  SCALING 


As  GFI  short  wavelength  experiments  were  not  done,  no  experimental 

information  is  available  on  wavelength  scaling.  Theoretical  estimates  of 

2 

hot-electron  temperature  tend  to  scale  with  the  combination  1^*  (1^  - 

intensity,  A  «  wavelength),  indicating  much  weaker  effects  at  short 
wavelength.  There  is  no  theoretical  information  available  on  fast- 
electron  current  scaling  other  than  our  calculations  of  Sec.  3.4,  but  at 
short  wavelength  the  temperature  Is  too  low  to  provide  significant 
electron  escape,  so  the  current  should  drop  approximately  as 


exP-[*/Amfp(Th)]  , 


(5.9) 


with  l  the  density  scaleheight,  and  with  the  mean  free  path  A^ 
2 

proportional  to  T,  ,  presumably 
n 


(li2)2/3  . 


(5.10) 


Electric  and  magnetic  fields  on  the  laser  pulse  timescale  should 


vary  In  roughly  the  same  way  as  the  currents,  while  fields  at  higher 
frequency  may  scale  differently  with  wavelength  because  their  source  Is 
fluctuations  In  the  space  charge  layer  and  the  current.  The  wavelength 
scaling  of  the  high-frequency  fluctuations  Is  unknown.  The  non-slnusoldal 
"virtual  cathode  oscillations"  of  the  space-charge  barrier  near  the  target 
give  a  radiation  field  (Eq.  3.121) 


E 

vc 


”h(0)T 


-1 

r 


/Tf(X) 


at  the  fundamental  frequency  (~l/2  GHz,  Eq.  <3.113>). 

This  Is  to  be  compared  with  the  'quasistatic'  (~0.1  GHz)  field  due 
to  the  space  charge  of  the  longer  plume  of  escaping  fast  electrons  (Eq. 
3.122),  which  should  scale  as 


V«~  *  V-V’  -  <VTh-x>  • 

With  1^  ~  1A  and  e$  ~  1  keV  as  observed,  this  lower-frequency  field  at  the 

D  detector  (r  ■  1  m)  would  be  of  order  0.1  V/cm  for  a  20  era  plume,  and 

2  1/3 

would  have  a  risetime  of  order  2-10  ns.  Assuming  e$  «  T.  «  (IX  )*,J 

m  h  L 

2  2 

from  Eq.  (3.9),  assuming  L  *  vr ,  i.e.  ,  L  *  t  e$  ,  and  assuming  current, 

IB 

2  3/2 

I^,  scaling  as  (l^X  )  as  observed,  the  'quasistatic'  plume  electric 


field  would  scale  as 


with  t  the  laser  pulselength,  whereas  the  radiation  field  from  the 
virtual-cathode  oscillations  appears  to  scale  as 


E 

vc 


r 


(V2)"2 


> 


f 

osc 


C^x2)3'4 


with  the  laser  pulse  risetime.  Based  on  the  observed  plume  currents 

(~  1/2  A)  it  would  appear  that  the  fundamental  (~  GHz)  field  from  vc 

oscillation  has  smaller  magnitude  (~  0.5  mV/cm  at  laser  intensities 
10  2 

~  10  W/cra  )  than  the  monopole  field  of  the  plume,  unless  the 
effective  plume  length  is  ~2  cm  or  smaller. 


5.3  PRESSURE  SCALING 


From  10"6  to  10“^  torr,  fast  electron  currents  should  scale  only 
slowly  with  pressure,  the  dependence  being  due  to  (a)  collisional  slowing 
of  electrons  in  the  gas,  and  more  importantly,  (b)  reduction  of  the  space- 
charge  barrier  by  gas  ionization  near  the  target.  Secondary  electron 
currents  neutralizing  in  part  the  fast  electron  current  scale  as  pressure 
times  fast  electron  current,  but  have  a  large-scale  oscillation  at  late 
time,  the  amplitude  of  these  oscillations  depending  on  the  laser 
pulselength  and  the  secondary  current,  and  the  period,  as  estimated  in  Eq. 
(3.105),  depending  on  pulselength  t,  fast-electron  current  1^,  and 
pressure  P  as 


t 

08C 


(IbTP> 


-1/2 


(5.11) 
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The  net  current,  seen  as  low-frequency  (t  ~  t)  magnetic  field  via 
the  field  magnetic  monitor  (FMM) ,  should  scale  as 

Xnet  *  V1  “  con8t*  P  Ibt2}  •  (5.12) 

because  the  secondary  current  is  proportional  to  the  number  of  secondary 
electrons  («  I^Pt)  times  the  E2  field  («  1^)  times  the  acceleration  times 
t,  at  low  pressure.  The  early-time  history  of  the  current  should  be 
unaffected  by  P  because  secondaries  will  not  have  had  time  to  form  or 
move. 

Above  10”-*  torr,  the  fast -electron  slowing  reduces  the  size  of  the 
current  plume  (see  Sec.  3.7)  and  the  secondaries  become  less  inertia- 
limited  and  more  collision-limited 

(J  «  E  Instead  of  J  *  E)  . 
s  s 

Experimental  results  for  pressure-dependence  of  the  FMM  current 
(proportional  to  Inet)  were  shown  in  Figure  4.17  and  are  subject  to  a 
major  uncertainty  in  the  pressure  scaling  because  of  an  order-or -magnitude 
disagreement  in  the  pressure  reading  from  two  different  diagnostics. 

At  higher  pressures,  of  order  0.1  torr,  where  the  collision 
frequency  of  secondary  electrons  becomes  comparable  with  their  plasma 
frequency  or  with  the  beam  timescale,  the  large  oscillation  of  the 
secondary  cloud  is  damped,  and  both  E  and  B  measurements  should  and  do 
confirm  this  (Figures  4.15,  4.17).  This  leads  to  the  observed  20%  per 
decade  of  pressure  reduction  in  the  field  strength. 
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APPENDIX  A 


HOT-ELECTRON  TRANSIT  TIMES  IN  A  SELF-CONSISTENT 
GROUNDED-DIPOLE  POTENTIAL 


1.  Transit  time  to  r. 


,(V 


-  /  zr  ■  / 


dr 


0 

*  v0 


u(r) 


nh(0) 


2U 


Ut  *  '  T~  "  “TT  "  "  *  •  Let  L  =  tn  [  j]  • 
n  h  n  q 


From  ♦  “  InC  -  2  In  cos(k  -  ax)  with  ie  »  cos  ^/C,  a  =  /C/2 ,  and  with 
£nC  =  -L,  we  get 


♦  ■  -  L  +  2  In  sec(K  -  ax)  , 


where  x  =  r/*D  and 


T  1/2 

(— r — ) 

4we  n,  (0) 
n 


Here,  again,  u  ”  ~  u  “  (~  T  )♦  . 

tn  m  n 


v2/(—  T  ) 

Let  X  ■  e  ^  m  h  so  that  -  (—  T.  )  tnx . 

u  ran 


Then 
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i 


J 


£ 

2 


BaaaoMBaMMa6iCMMM!iflMM«iii4aai«ffliicaaaa«^^ 


y  Vq  -  u  -  (-  Th)  /*nx  -  +  tn  sec  (k 


-  (~T.)1/2  ^ ln[ R  sec2(tc  -  o> 

n  n 


where  R  -  xC  >  1  for  transmitted  particles. 
Thus 


T<V  - 


tc/a 


dx 


0  /ln[R  sec^(ic  - 


since  x(r  )  -  </a  . 

q 


,  x  .-1/2  -1  ( 
r(r  )  -  2  u  J 

q  p 


x/a 


dx 


^  )J ln{R  sec2(x  -  a; 


where  <o  ■  (4*e^n  (0)/m)^^. 
P  *» 


To  do  the  Integral 


i  -  r 


x/a 


dx 


°  An[R  8ec2(ie  -  ax)  ] 


let 


\/ln(R  sec2(K  -«*x)) 


so  that 


Upper  Llalt:  see  (0)  -  l  ,  w  «  /inR 


Lower  Llalc:  (*e  -  cos_1/N  sec2*  -  1/N).  w  -  /tn(R/N). 


x  ri/2  /I^r7nT  dw 

“  /S3  JV~ 

*  e  -  R 


It  Is  this  Integral  we  now  evaluate. 


Let  w  «yln[R  sec  (*  -  ox)],  R  >  1,  as  before. 


Note  that  Min.  w  -  / InR ,  at  x  ■  x/a, 


2  2 

Let  wQ  s  /InR.  Expand  ew  about  w*  «  InR. 


Let  w  =  u  and  u^  -  in  R. 


urt  +  Au  u_  Au  f. 

u  0  0 

e-e  -  e  e  -  Re 


R [ 1  +  (u  -  Uq)  +  j  (u  -  uQ)2  +  £  [u  -  uQ)3  +  ...]  . 


eW  -  R  =  Jeu  -  R 


R1/2(u  -  urt)1/2  [1  +  i  (u  -  uQ)  +  |  (u  -  uft)2  +  ...J1/2 


R 


1/2  _ dw 


-  r-~ - m  [  1  -  |  (w2  -  InR)  +  ^  (w2  -  tnR)2  +  . 

(t»  -  tnR)1' 


f72 

*  e  -  R 

on  expanding  the  [  |*/2  square  root. 


•  •  •]  > 


But 


l 


tn(R/n)  - 


/tnR  + 


tn( 1/N) 


+  L 


(L  >  1) 


So 


Rl/2/ 


/TOTnT 


dw 


dw 


TnR 


/  2  /I  2 

*ew  -  R  /w  -  wQ 


-  f  /  (w2  - 


2\ 1/2 
wQJ  dw 


+  ^/  (w2  -  w0)3/2dw 


f~2 

with  lower  limits  wq  and  upper  limits  /Wq  +  L. 
First, 


/ 


dw 


cosh 


Ju 


■F5 


’“2  -  w0 

2  ,  .  ,  f  Vt  Th^ 

But  wQ  -  tnR  -  tn(xc)  “  tn[e 


N] 


(>o) 


-  L  . 


2  2  2 

When  vrt  -  l(—  T, )  =  v  ,  particles  are  Just  barely 

u  ran  € 

transmitted  over  the  barrier.  So  in  terras  of  the 


escape  velocity  v£ , 


wo  TJ 

Oran 


voL(Z  v 


V-  -  V 
0  € 


2 

2  r  v°  1 1 

v 

e 


.  2  ve  2.  .  a  -  rv°  n 

l.e. ,  w  - —  -  v  A,  where  A  =  {— -  1 )  . 
u  l  e 
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Thus 


U  m  LA  • 
0 


rw0  +  L 


rt — s' 

v  w  -  w0 


-1  /TT?. 


Call  this  I, 


Next, 


*  ~tI 


<v0  +  L 


rw  -  Wg  dw 


-  -  |  {✓  1  +  A  -  A  cosh-1  J 1  +  1}  . 


Here  L  - 


tn(l/N)  and  -  LA,  so  /w^  +  L  /L  ■ 


1  +  A  . 


Likewise, 


,  1  r  "  (2  2>3/2, 

+  J  lw  -  wQJ  dw 


V  {/rTir  -  7  LyirrT  +  7 


Combining  Ij  terms, 


ll  +  *2  +  X3 


. -l  /.  \i  r.  ^  i  r1^  i  rLA^2 
cosh  /l  +  T  •  [  l  +  y  ( j-J  +  -jj  (— )  + 


+ 


] 


-  ^  ♦  is  -  X>2 

(Example:  For  L  -  4  and  A  -  1,  Ij  -  .881,  I2  -  -.266,  I3  -  .0256, 
X1  +  X2  +  X3  "  *64)* 


Then 

T(rq)  -  2-1/2a-1uT1[I1  +  I2  +  I3  +  _ ]. 

But 

a  -  /C72  so  r(rq)  -  (Cu»2)_l/2  [I }  +  i2  +  13  +  . . .  ]  , 


inversely  proportional  to  the  1/2  power  of  the  hot  density  at  rq. 


2.  Transit  time  from  r^  out  to  r  »  r^: 

Next  we  use  a  model  for  the  potential  at  r  >  rq  to  calculate  the  rest 
of  the  transit  time  out  to  r  »  Tq, 


T 


2 


(x  -  t(r  )  +  r_)  . 

q  Z 


The  potential  model  must  have  the  value  U(Tq)  at  r  *  Tq,  matching  the 

2 

value  of  the  inner  solution  (r  <  rq),  ar,d  n^st  go  to  zero  as  1/r  as 
r  ♦  •  in  order  to  be  dipole-like.  We  take 


U  -  U(rq) [ 1  +  (r  -  rq)2/r2] 
for  r  >  r  .  Let  s  -  (r/r  )  -  1.  Then  (r  -  r  )2/r2  -  s2. 

q  q  q  q 

Again  let  u  ■  2U/m,  u(s_0)  ■  v2,  and  A  ■  (vq/v2)  “  1  so  that 
V0  ‘  V£(1  +  A)* 
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m 


Then 


r  X-l 

t2  *  ^  / 
z  ve  0 


i  +  r  d< 


r(l  +  A)(l  ♦  s‘)  -  1 


r  X-l 

-  -*/ 


1  +  s‘ 

A  .  2 

1  ♦  A  * 


where  X  Is  s  +  1  evaluated  at  r  -  detector  radius,  i.e. ,  X  ■  r,  _/ r  . 

aec  (| 


O  I 


Since  the  integral  is  not  known  in  closed  form,  we  break  the 

1  X-l 

integration  region  further  into  /  ds  and  /  ds.  In  the  former,  we  expand 

0  1 

in  X  .  In  the  latter,  we  expand  in  1/X^. 

For  s  <  is 


1  ♦  s 


1  +  A, 1/2 


■  (-4-=o 


i  1  .  * 

1  +  3 

i  ♦  (-Lfi) 


-  (  )WIU  -  j  +  2(y2  +  y3>W  +  *)  -  ...)  , 


where  y  -  sV2A. 

This  seeos  to  converge  only  on  x  <  so  to  have  it  applicable  for 

x  <  1,  we  must  have  S7X  >  1,  i.e. ,  A  >  1/2. 


For  the  s  >  1  integration,  since  the  integrand  is  nearly  unity,  we 
use  the  fact  that 


where  6  =  A/(l  +  A),  i.e.  ,  subtract  off  the  undelayed  l/v^  transit  time 
from  the  integration,  adding  it  back  outside  the  integral.  Expanding  the 
numerator  (only)  of  the  new  Integrand  in  1/s^: 


Thus  we  have 


'"geas  [  jn? 

1  V  5  +  8 


-  l]  *  ?  I  (-D 


k  (1/2  -  k)! 


k-1 


k! 


(1  - 


6k)  / 


c,. 


(y 


This  can  be  written  as 


(26)  1/2{  /T7I  A(6)  +  [Lz±  +  6B(6)]tn  (^IZZjL^)} 

✓  1+6  -  /6 


where 


A(6 )  = 


y  y  (-i)j^1  k){  61-J 

k-2  (  U  k! (k-j) 


j-1 


(1  -  6^)  IT  (l  - 


2(k 


and 


•  ..  (7  -  k)!  .  k-1 

B(6 )  =  l  (-l)2k+1  -  (6-k  -  1)  n  (1-i-) 


k=l 


j-1 


This  Is  derived  using 


i  .  I  r/1  +  « 

1,  A  l 


6  L  k  -  1 


^  “  2(k  -  D^k-J 


and 


1  .  //T+T  +  /T, 


T(2r  to  r) 

q 


Is.  JL 
v0  /26 


1  +  j  A(6)  + 


/ 1  +  5  -  /6 


i 
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APPENDIX  B 


HOT-ELECTRON  TRANSIT  TIMES  IN  A  MONOPOLE 
SPACE-CHARGE  POTENTIAL  HELL 


For  a  potential  as  in  Flg.Bl  Hue  Co  a  charge  q  inside  radius  rq,  one 


has 


2  tc  jc 

rQ  +  r  -  with  vf 

q 


^•c'/r  the  escape  velocity; 

r»g.  bi. 


1  •  c«  • 


A-*-*  • 


2  2 

Let  r/r  -  x,  and  again  let  A  •  Vq^  ”  1  •  fraction  of  energy  above 


e  .  Then 
esc 


r  ,  x 

fc  *  ■  r,/ 


dx 


i  n. — rr — rr 

/v0  -  ve  (1 


-  A-  {^n  »Yx7  -  '/TTI  ♦  A"1/2tn  1  |}  . 

v  ♦  a/71 


Since  vo  "  ve  "  veA’  w*  hava  v0  "  ve(I  +  and  thu*  ^  ^  — ■ 


v0  ve  /TTT 


For  large  A  -  1  (vg  >  ve ,  early  arrivals)  at  x  »  1: 


/'  |x  -  -V  {i1/2<*  -  1)  ♦  4-1/2tn  SS-  1 

re  A/  x 
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with 


M  =  1  +  /T  +  A 


=r 


1  +  i 1  +  (Ax)"1 


~  i 


and 


/ 1  +  A 


So 


J  r 
r 

q 


Is.  I  1  _  .  J 1  +  s'1 

it  I*  -  24  +nx  +  — 


r  r  1  fnx  N  -  1 
v  * 1  “  2A  x  2Ax 


2A 


1  +  A 


where  N  =  / 1  +  A  1  ~  1  +  (1/2A). 

^1 

Here  the  lnx/2Ax  Is  small  for  large  x  (and  If  A  ,  the  N  -  1  term  is  even 
smaller).  Note  that  and  r/vQ  is  the  unimpeded  x. 

Thus, 


and 


T  “  ~  { 1  -  2^7  for  large  A  and  large  x, 

T  ’  ^  7T7T  t„  |/Ax  +  /T  +  A*|} 

for  small  A  and  large  x. 


i.e. ,  for  electrons  hear  the  escape  energy. 
From  this. 


dx  _  x_  2r_  9JF 
dV0  ~  v0  v2  3A  ’ 
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In  the  first  case 


or 

F  -  '  {--X(y--— ^  *n|^x  +  ’'T^nnri} 


In  the  second  case 


In  the  first  case: 


3F  lax  1 

jk  “  tzt  ~j 

A 


for  large  A 


so 


dr 

dv 


0 


-  2_  +  £_  r  &nx  l  +  a-i  l 

t«  t#  ^  ^ J  y 


V0  V0  *  A 


r  rlnx  1  +  A  ,  ,  tnxi 

l-= - t-  -  1  + 


~T  - 7? 

v0  A 


r  r r 3  .  l\  inx  _  ,i 

T  Ll7  F  “57 


-1  «1 


In  the  second  case: 


1/2  j/Ax(l  +  Ax)  -  £n  |/Ax  +  /l  +  Ax 


{ 


(Ax) 


!7T 


:112  {^LZ  -  1  In  | y1/Z  Wl  +  y|  }  (y  =  Ax)  . 

y  yJ/ 


So 


l_(r/„,/2>  - 


,  3  I  -  |  1/2  ,-r— - 

T  “T7T  £n  ly  +  /I  +  y 


T77  f;‘«  ^ + 


which  approaches  2y“^^  for  large  y. 


For  small  A  but  x  »  1,  with  Ax  ~  1  (Middle  arrivals),  we  have 


JT  dr 
r  v 

q 


1  +  Ax  -  (l  +  i)  +  A_1/2£n  I  - U-LA -  1} 

1  /A  /l  +  Ax  +  A/x  1 


JL  ltin  /Y~*~Sx  -  (I  ♦  *)  *  *-*«».  I  ■  1  ^  *‘/z  1} 

0  '  /I  +  Ax  +  /AT  1 


Yf 2 —  {*^Ax(l  +  Ax)  -  £n  |/Ax  +  /l  +  Ax|  +  £n(l  +  /A)} 


vqA  Ax 


When  Ax  is  fairly  large,  this  is 


r  dr  rq  x  rl  , 

Jr? - T77  "  v - T7T  *  neglecting  the  £n  term. 

V0  L  '  0  L  ' 


APPENDIX  C 

EPPECT  OF  LOCAL  MAGNETIC  FIELDS  ON  ELECTRON  PLUME  MOTION 

It  is  instructive  to  look  at  the  typical  electron  magnetization  and 
E  *  B  drift  in  the  magnetic  fields  due  to  Vn/n  x  7T  generation  as  well  as 
due  to  the  fast-electron  current.  Whereas  the  beam  self  field  is  of  order 


B(G)  -  0.2  Ib(A)/a(cm)  , 


(C.l) 


i.e. ,  a  few  Gauss  at  most,  giving  a  gyroradius  of  £  10  cm  for  1  keV 
electrons,  by  contrast  the  more  localized  fields  due  to  Vn  x  VT  have  much 
higher  peak  values  in  the  ablation  layer.  From 


3tB  -  (c2/4iro)V2B  -  V  x  x  B) 


Vn 

S.  — 1  x  VT 
e  n 

e 


(C.2) 


one  has  near  the  edge  of  the  heated  spot  (where  Vn  and  VT  are  most  nearly 
perpendicular)  a  B-field  generation  rate  of  order 


B 

x 


X  VT 


lo8Vh  A 


(C.3) 


with  hn  and  hT  the  density  and  temperature  scalelengths  (cm)  and  in 
G/s.  Neglecting  the  saturation  of  this  field  build-up  for  the  moment,  one 
estimates  at  early  times 


B(G) 


t(n8) 


200  •  T  (^£5cm)(i00li“)  F(*)  , 

ev  n 

i  n 


(C.4) 
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where  F(r)  has  maximum  value  unity  and  decays  rapidly  with  distance  from 
the  edge  of  the  laser  heating  spot.  For  T  ~  10  eV,  this  would  give 

B  ~  20  kG  •  F(r)  in  10  ns  (C.5) 

if  the  build-up  of  B  were  not  saturated  earlier.  The  local  fast-electron 
gyroradlus  in  such  a  field  is 


rg(cm)  -  75.3  [Th<keV) ]  w  VB(G)  , 


(C.6) 


which  could  decrease  to  r  ~  30  um;  and  the  gyrofrequency 

ty 


u  (s'1)  -  1.76  X  1010  B(kG) 

ce 


(C.7) 


could  give  gyroperiods  £  2  *  10  s,  i.e. ,  perhaps  short  enough  for  the 


hot  electrons  to  be  well  magnetized  until  they  drift  out  of  the  high-field 


region. 


The  cE/B  drift  energy  of  an  electron  is  given  by 


jmv^(eV)  -  2.7  [ 


E(kV/cra)i2 

I'M'"'1 


(C.8) 


when  this  is  less  than  T^> 


If  E  is  derived  from  a  potential  drop  of  order  1  keV  over  a  sheath 


scalelength  of  order  10”^  cm,  the  square  bracket  can  be  of  order  unity 


where  E  and  B  are  largest.  Drifting  at  this  rate,  a  plume  electron 


escapes  the  high-field  region,  of  size  t  in  a  time 


k. 


I 


s 


y 

y 


pi* 


I 


| 

Eft 


'*’»  .‘'V  ’V  "V*  » '  ji**  -Nyv*  *^**j',*>. 


t  -  ~  10  ns  •  t(c»)  •  (C.9) 

When  £  Is  microscopic,  e.g. ,  100  urn,  this  Is  often  negligible;  when  It  Is 
macroscopic,  e.g.,  1  cm,  electron  densities  and  transit  times  are  affected 
on  the  experimental  timescale. 

But  as  the  magnetic  field  Is  generated  at  the  edge  of  the  heating 
spot,  It  diffuses  and  propagates  away,  so  that  during  the  pulse  the  local 
peak  value  of  B  usually  Is  limited  to  a  value  below  that  indicated  in  Eq. 
(C.4). 

If  the  copper  target  conductivity  Is  ot  and  its  thickness  £t,  and  if 
the  ablating  plasma  has  conductivity  cr  and  thickness  £  ,  one  can  estimate 
for  the  limiting  B, 


4*T 


c  1 


1 


max 


V'  [„\2  +  o  t2  ♦ .  t2 

a  a  a  a  t  t 


-1 


(C.10) 


where  6r  Is  the  scalelength  for  nonuniformity  across  the  spot;  this  limit 
Is  reached  in  a  time 


<  ~r 


1 


1 


-1 


o  £2 
a  a 


2  2 1 
o  £  +  o  £ 

a  a  t  t 


(C.ll) 


For  copper,  a  ~  4  x  10^  e.s.u.  ,  and  the  target  was  thick 
-  £  ~  0.63  cm.  For  a  10  eV  plasma  blow-off  with  l equal  to  the  density 


gradient  scalehelght,  of  order  100  inn,  one  has  a  Spitzer  conductivity 


2  9  2  2 

giving  cl  *  9  *  10  c«  s.  Thus  o  t  can  be  neglected  compared  with 

as  a  a 


2  O 

a  1  ,  and  the  square  bracket  in  Eq.  (C.ll)  is  dominated  by  o  is 

t  t  ad 


B  <G>  ~  0.25  (T/10  eV)(a  /  9  x  1013s_1 )(i  /100  pm)(.05  cm/6r) 

IBaX  a  8L 


tg  ~  0.13  ns  . 


(C.13) 


So  B  may  be  limited  early  in  the  pulse  to  peak  values  which  increase 
quasi-statically  as  T5/2^~  250  G,  giving  r  ~  .2  cm  for  the  hot  electrons 


(gyroperiod  ~  1  .A  ns),  and  in  this  case  the  fast  electrons  are  not 


strongly  magnetized  and  probably  leave  the  strong-field  region  in  a 


portion  of  a  gyroperiod  rather  than  by  drift  motion. 
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Conversion  factors  for  U.S.  Customary  to  metric  (SI)  units  of  measurement 


MULTIPLY  - ►  BY  - *•  TO  GET 

TO  GET  ^ -  BY  ^ -  DIVIDE 


angstrom 

1. 000  000  X  E  -10 

meters  (m) 

atmosphere  (normal) 

1.013  25  X  E  42 

kilo  pascal  (kPa) 

bar 

1.000  000  X  E  42 

kilo  pascal  (kPa) 

bam 

1.000  000  X  E  -28 

2  2 
meter  (m  ) 

British  thermal  unit  (thermochemical) 

1.  054  350  X  E  43 

joule  (J) 

calorie  (thermochemical) 

4.  184  000 

joule  (J) 

cal  (thermochemical) /cm 

4. 184  000  X  E  -2 

mega  joule/m2  (MJ/m2) 

curie 

3.  700  000  X  E  4l 

•gig a  becquerel  (GBq) 

degree  (angle) 

1.  745  329  X  E  -2 

radian  (rad) 

degree  Fahrenheit 

tg  =  (t*f  4  459.  87) /l.  8 

degree  kelvtn  (K) 

electron  volt 

1.  802  19  X  E  -19 

joule  (J) 

ei* 

1.000  000  X  E  -7 

joule  (J) 

erg/second 

1. 000  000  X  E  -7 

watt  (W) 

foot 

3.  048  000  X  E  -1 

meter  (m) 

foot-pound  -force 

1.  355  818 

joule  (J) 

gallon  (U.  S.  liquid) 

3.  785  412  X  E  -3 

3  3 

meter  (m  ) 

inch 

2.  540  000  X  E  -2 

meter  (m) 

jerk 

1.000  000  X  E  49 

joule  (J) 

joule  Ailog  ram  (J/kg)  (radiatton  dose 
absorbed) 

1. 000  000 

Gray  (Gy) 

kilotons 

4.183 

terajoules 

kip  (1000  Ibf) 

4.448  222  X  E  43 

newton  (N) 

kip/lnch2  (ksi) 

6.  894  757  X  E  43 

kilo  pascal  (kPa) 

ktap 

newton  -second  /m 

1.  000  000  X  E  42 

(N-s/m2) 

micron 

1  000  000  X  E  -8 

meter  (m) 

mil 

2  .  540  000  X  E  -5 

meter  (m) 

mile  (international) 

1.  609  344  X  E  43 

meter  (m) 

ounce 

2. 834  952  X  E  -2 

kilogram  (kg) 

pound-force  (lbs  avoirdupois) 

4.  448  222 

newton  (N) 

pound -force  inch 

1.  129  848  X  E  -1 

newton-meter  (N-m) 

pound  -force /inch 

1.  751  268  X  E  42 

newton/meter  (N/m) 

pound -force  Aoot 2 

4. 788  026  X  E  -2 

kilo  pascal  (kPa) 

pound-force/inch2  (psi) 

6.  894  757 

kilo  pascal  (kPa) 

pound-mass  (Ibm  avoirdupois) 

4.  535  924  X  E  -1 

kilogram  (kg) 

pound -mass -foot2  (moment  of  inertia) 

2 

kilogram-meter 

4.  214  011  X  E  -2 

0«-m2) 

pound -mass  Aoot2 

kilogram/meter 

1.  601  846  X  E  4l 

(kg/m3) 

rad  (radiation  dose  absorbed) 

1.  000  000  X  E  -2 

••Gray  (Gy) 

roentgen 

2.  579  760  X  E  -4 

coulomb  Ailog  ram 

(C/kg) 

shake 

1 . 000  000  X  E  -8 

second  (s) 

slug 

1.459  390  X  E  4l 

kilogram  (Ig) 

torr  (mm  Hg,  0*  C) 

1. 333  22  X  E  -1 

kilo  pascal  (kPa) 

•The  becquerel  (Bq)  is  the  St  unit  of  radioactivity;  1  Bq  *  l  event/a. 

••The  Gray  (Gy)  is  the  SI  unit  of  absorbed  radiation. 

A  more  complete  listing  of  conversions  may  be  found  tn  "Metric  Practice  Guide  E  380-74  •• 
American  Society  for  Testing  and  Materials. 
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